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Abstract: We continue our study on the gravity duals for strongly coupled large 
p i' N QCD with fundamental flavors both at zero and non-zero temperatures. The 

^ ■ gravity dual at zero temperature captures the logarithmic runnings of the coupling 

constants at far IR and the almost conformal, albeit strongly coupled, behavior at 

■ the UV. The full UV completion of gauge theory is accomplished in the gravity side 

■ by attaching an AdS cap to the IR geometry described in our previous work. At- 
' taching such an AdS cap is highly non-trivial because it amounts to finding the right 

Q ■ interpolating geometry and sources that take us from a gravity solution with non- 

■ zero three-form fluxes to another one that has almost vanishing three-form fluxes. 
In this paper we give a concrete realisation of such a scenario, completing the pro- 
gram advocated in our earlier paper. One of the main advantage of having such a 
background, in addition to providing a dual description of the required gauge theory, 

^ ■ is the absence of Landau poles and consequently the UV divergences of the Wilson 

■ loops. The potential for the heaviest fundamental quark anti-quark pairs, which are 
like the heavy quarkonium states in realistic QCD, can be computed and their linear 
behavior at large separations and zero temperature could be demonstrated. At small 
separations the expected Coulombic behavior appears to dominate. On the other 
hand, at non-zero temperatures interesting properties like heavy quarkonium type 
suppressions and melting are shown to emerge from our gravity dual. We provide 
some discussions of the melting temperature and compare our results with the Char- 
monium spectrum and lattice simulations. We argue that, in spite of the large 
nature of our construction, certain model-independent predictions can be made. 
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1. Introduction 

The study of strongly interacting matter under extreme conditions of temperature 
and/or density is one of the most fascinating areas of contemporary subatomic 
physics. This program aims to explore the many facets of the bulk behavior of 
Quantum ChromoDynamics (QCD): The theory of the strong interaction. It seeks 
to map out the different phases allowed by QCD and the nature of the possible phase 
transitions connecting them or in short, to elucidate the QCD phase diagram. In 
this context, the existence of an exotic phase of QCD, a quark-gluon plasma, has 
been a prediction of lattice QCD whose details have continuously being refined over 
the years [|l|. On the experimental front, several observables have been put forward 
as signature of the quark gluon plasma. These include electromagnetic radiation [0] , 
the quenching of energetic QCD jets 0, and the dissolution (with increasing collision 
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centrality and energy) of heavy quark bound states according to the seminal sugges- 
tion in Ref. 0]. The Relativistic Heavy Ion Colhder (RHIC), now at the end of its 
first decade of operation, has uncovered an intriguing set of phenomena suggestive 
of new physics. One of these is the observation of strong hydrodynamic flow effects, 
highly suggestive of a "strongly coupled" quark gluon plasma p. 

The fate of quarkonium is being analyzed at RHIC, as it was at the SPS before. 
A surprising fact to emerge of these studies is that the suppression of the cc ground 
state - the J/ip - at RHIC is entirely comparable to that at the SPS, in spite of 
the much larger energy densities being reached at the first facility. This triggered 
many analyses with scenarios where the enhanced dissociation at RHIC was roughly 
compensated by an extra formation owing, for example, to quark-antiquark coales- 
cence near hadronization 0. Related investigations are concerned by the fate of the 
quarkonium spectral density above Tc, the deconfinement temperature 0. It is fair 
to write that the study of quarkonium imbedded in a finite-temperature strongly 
interacting medium is a flourishing industry: The modifications of its spectral profile 
can be related to in-medium effects. A related topic of investigation on the lattice 
consists of calculating the quark-antiquark potential as a function of temperature [§]. 
These calculations show a Coulomb potential at zero temperature, with an added 
linear part that slowly disappears as T is raised, leading eventually to the unbinding 
of quarkonia bound states. Our goal in this work is to approach this softening of the 
potential from a different point of view. 

In parallel with the studies described in the previous paragraph, the physics of 
hot and dense strongly interacting matter, and thus that of the quark-gluon plasma, 
has recently benefited from the use of a new set of techniques, germane to string 
theory. The gauge-string duality can indeed provide a sophisticated toolbox with 
which to treat strongly-coupled, strongly interacting systems [|], |10], |ll|]. Our purpose 
here is to bring closer the more traditional investigations in QCD with those pursued 
in string theory. In gauge-string duality, a finite-temperature medium is dual to 
a black hole. Even though in a large number of applications the associated field 
theory is conformal, we use a framework which is "QCD-like". More specifically, we 
construct the dual gravity of thermal field theory which becomes almost conformal 
in the UV but has logarithmic running of coupling in the IR with matter in the 
fundamental representation. Without being explicitly QCD, this string theory will 
provide some of the features associated with large N Quantum Chromodynamics, 
and its study may shed more light on the behavior of strongly coupled, strongly 
interacting matter at finite temperature. 

Our paper is organized as follows: In the next section we define the geometry in 
which our solution will exist. The description of the full geometry is subtle, so we 
will divide the geometry in three regions. The far IR geometry will be described in 
sec. 2.1, and the far UV geometry will be described in sec. 2.3. These two geometries 
are connected by an interpolating geometry that we will describe in details in sec. 
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2.2. Once we have the full geometry, we compute the heavy quark potential from 
the Nambu-Goto action first for zero, and then for finite temperature in sec. 3. 
In this section we will also provide a generic argument for confinement both for 
zero and non-zero temperatures. Although most of our analysis in this paper will 
be done analytically, we will do some detailed numerical analysis to study regimes 
that are difficult to access analytically. We will show that the numerical analysis 
fits consistently with the expected behavior of the heavy quarkonium states in this 
theory. Finally, we summarize and conclude. 



2. Construction of the Geometry 



Following the development in |0, in [|T3| it was shown that a geometry where the 
dual thermal field theory was almost conformal in the UV, and had a logarithmic 
running of the coupling in the IR existed. The gauge theory studied in |TB[ had a 



dual weakly coupled gravity description at zero temperature in terms of a warped 
deformed conifold with seven branes and fiuxes. The gauge theory in turn is strongly 
coupled with a smooth RG fiow but no well defined colors at any given scale. When 
the gauge theory is weakly coupled, the description can be presented in terms of 
cascades of Seiberg dualities that slows down quite a bit when one approaches the 
far IR because of the presence of fundamental fiavors. There is no supergravity dual 
description available for this case, and the cascade is only captured by the full string 
theory on the relevant geometry. 

Once a non-zero temperature is switched on, the strongly coupled gauge theory 
description is given by a dual supergravity solution on a resolved warped deformed 
conifold with seven branes and fiuxes [jl3l. The resolution factor is directly related 



to the temperature because in the presence of a black hole a consistent solution of 
the system can only be achieved by introducing a non-zero resolution factor for the 
two cycle. In a Klebanov-Tseytlin type geometry, this resolution factor would in fact 
remove the naked singularity. 

One of the other key ingredient of the solution presented in is the far UV 



picture. In all the previous known attempts to this problem, the dual supergravity 
solution was always afflicted by the presence of Landau Poles. Such problems arose 
because of the behavior of the axio-dilaton, that typically blows up due to their 
logarithmic behavior. What we pointed out in |jl3[ is that the logarithmic behavior. 



which is so ubiquitous in these constructions, appears because we are studying the 
theory near any one of the seven branes. In the full F-theory picture the large r 
behavior is perfectly finite, and in fact also has a good description in terms of the 
metric too. The behavior of the warp factor for large r is given by: 

" = (2-1) 



3 



with r(a) = r^^ ~ and e(^a) is a small positive number that is a function of ggNf, QsM 



and QsN (see eq. (3.36) of The sign difference from \Ti\ is just a matter of con- 
vention). The log r appearing in the warp factor doesn't create much of a problem 
at UV: the theory is perfectly holographically renormalisable, and any fluctuations 
of the background are under control. The fact that we can have well defined and 
renormalisable interactions in this background in the presence of fundamental fla- 
vors was shown, we believe for the first time, in ||13[ (see ||14| for renormalisability 



argument without fundamental flavors). For the present purpose, we want to ask a 
slightly different question here, namely: can we construct a dual supergravity back- 
ground that allows logarithmic RG flow in the IR but has a vanishing beta function 
at far UV? From our discussion of the UV caps in it is clear what we should 



be looking for: we need a gravitational background that resembles OKS geometry 
for small r, but has a UV cap given by an asymptotic AdS geometry. To extend 
this configuration to high temperature, we need OKS-BH geometry^ at small r, and 
asymptotic AdS-Schwarchild geometry at large r. 

Such a geometry looks complicated, so we may want to ask whether we can switch 
off the three form fluxes and still have a dual description with running couplings. If 
this were possible then the analysis could be made much more simpler. It turns out 
however that such a simplification cannot occur in our set-up. To elucidate the last 
point, let us give a brief discussion. 

The RG runnings of the two gauge groups in this theory are determined by the 
following dual maps in terms of the bulk axio-dilaton r and NS potential B: 



47r2 47r2 



H TT- = TT Im r 



9l gl 

47r^ 47r^ Im r 



Im T f 

/ B-nlmr (mod 27r) (2.2) 

27ra' Js^ 



gf gi 2'na' J 32 

Once we switch off B the two couplings would be the same and would induce a 
Shifman-Vainstein /3-function of the form: 

r\ Q 2 

3N-2N{l-jA,B)-Nf{l-jg) (2.3) 



d log A g^^ 

where A is the energy scale that is related to the radial coordinate r in the gravity side, 
and Ta,b and 7^ are the anomalous dimensions of bi-fundamental and fundamental 
fields respectively. With such a picture of the flow, we might think that the F-theory 
completion might be to simply add sufficient number of seven branes parallel to the 
spacetime directions and wrapping the two internal two-spheres (so that they are 
points in the (r, ■?/') plane). This simple picture would unfortunately be inconsistent 



^For more details on the construction of OKS-BH (Ouyang-Klebanov-Strassler-Black-Hole) ge- 
ometry, see sections 3.1 and 3.3. 
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with the underlying cascading dynamics as could be seen from a T-dual framework 
0, 0, and therefore would be incapable of showing certain important behavior 
expected from this model. 



To understand the problem, observe that in the T-dual picture d la [|I6|, the 
D3-D7-conifold geometry is mapped to a configuration of D4-D6 and intersecting 
NS5 branes. The NS5 branes, that are T-dual to the conifold, are along x^^^^^^ and 
^012389 cJix-ections and have N D4 branes (T-dual of N D3 branes) between them. 
Everytime we cross the NS5 branes we expect extra D4 branes to appear because 
of the D6 branes. This is however only possible if the D6 branes are along 
which in turn would imply that in the brane side the D7 branes have to be along the 
radial direction. Additionally, motion of the NS5 brane would imply a, B^s held in 
the brane side that is not a constant but has at least a log r dependence along the 
radial direction. This means that Hj^s is non-zero, and we need to switch on H^r to 
satisfy the equations of motion, bringing us back to the model originally advocated 
0! 



m 



The discussion above should convince the readers that there aren't much avenue 
to simplify the original proposal of [|13] . The original model proposed in ||13| is struc- 



turally complicated, but is possibly the simplest in realising some of the properties 
of IR large QCD. A model simpler than this would be deviod of any interesting 
physics. 

Once this is settled, we want to see how to construct the kind of geometry 
that we mentioned above. Our requirement is to impose confinement at far IR and 
vanishing beta function at far UV. Since the original model studied in doesn't 
quite have the right large r behavior because the warp factor therein goes as we 



need to add appropriate UV cap. However before we actually go about constructing 
the background, let us clarify how addition of UV caps in general can change IR 
geometries. One thing should of course be clear, the far IR geometries cannot change 
by the addition of UV caps. This is because the UV caps corresponding to adding 
non-trivial irrelevant operators in the dual gauge theory^. These operators keep far 
IR physics completely unchanged, but physics at not-so-small energies may change a 
bit. So the question is how are these changes registered in our analysis? Additionally 
we may also want to ask how entropy of our gauge theories affected by the addition 
of UV caps? 

Both the above questions may be answered if we could figure out how the UV 
caps affect the energy momentum tensors of our gauge theories. The generic form 



of the energy- momentum tensor that we derived in our earlier paper [13] can be 



^Assuming of course that the relevant operators were responsible for creating the cascading 
dynamics from a given UV completion in the first place! 
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reproduced as: 



rpmm 
medium+quark 



+ (i^H + '1 + E ^S)^"'^-^"'^"^ + (^CTe--^-) (2.4) 



j=0 



where and -ft'i™" depend on the full background geometry via eq. (3.124) in 

^(i) 



[0 with s™'^' being the Fourier coefficients. The other terms namely and A/u- 



together specify the boundary theory for a specific UV completion [ [L3|] . 

Now its easy to see how the UV caps would change our results. Once we add 
a UV cap the local region — ai < r < + a2 near the junction^ at r = 
changes, with {ai,a2) being some appropriate neighborhood around r^. This means 



that Cl"",^™"' and i?™" etc. in eq (3.124) of would change. These changes can 
be registered as 



^iT ^ = ^iT + ((^^ua) -^^"a))e"^['"'('^']^°'' + C(e-^'^"^) (2.5) 

where the last terms in both the above equations appear from additional UV degrees 
of freedom, C]^^^ etc are the relevant a-th components of (7|"" etc, and Af^s is the 
effective number of degrees of freedom at the cutoff. 

Once we know these changes, its not too difficult to figure out the changes in 
the entropies due to the addition of UV caps. All we need are the RHS of eq (3.220) 
T3| using the results from (^75|) and taking care of the boundary temperatures 



from the changes in the warp factors'^. Using (p.5|) the result can be written as: 



s \r 2h(T) dT J 



+ 





\a\ |q! l^'^nn 


5(5)[/3]^ 








"o-(mn)~(4)[/3] ^(mn) ( ~S^m\ 


+ C(e--^"-) 



However physics that are only sensitive to far IR dynamics of our theory will not be 
affected by the addition of UV caps. On the other hand in all cases, far IR or not, 
none of our results could depend on the cut-off Tc- The results are only sensitive to 
the changes in IR geometries (via (p.5|) ) and the UV degrees of freedom (via e"-'-^""). 



^Clearly ai << Tc because the far IR geometry should remain completely unaltered. 

^There may be interesting cases where the changes in the energy-momentum tensors are com- 
pensated by the changes in the boundary temperatures. In such cases the entropies may remain 
unchanged. Here we will not consider such cases. 
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From the above discussions we see how IR geometries could be affected by the 
addition of UV caps. This then tells us that we cannot simply add an AdS geometry 
at r = Tc- The vanising beta function at UV could be realised by an asymptotic AdS 
geometry, a geometry whose warp factor behave as only asymptotically. In other 
words we require: 



^4 



.j,k=0 



bjk{ip,Oi,(j)i)\og''r 



for large r 
for small r 



(2.7) 



where {6i,(j)i,ilj) are the coordinates of the internal space. Observe also that we are 
now identifying the small r behavior of the warp factor to the relation (|2.1| ) given 
above. The precise connection will be spelled out in details below. 

Let us now make this a bit more precise. We require a gauge theory with confining 
IR dynamics and almost free UV dynamics at zero temperature, and then we want 
to study this theory at a temperature higher than the deconfining temperature, as 
mentioned before. Our dual gravity background that could in principle reproduce 
the gauge theory dynamics couldn't be the pure OKS (or OKS-BH) background of 
1^. We need an appropriate UV cap. Again, as we mentioned earlier, the UV cap 



should be asymptotically AdS. The warp factor should have the form ( PT7| ) at UV 
and IR, so we need an interpolating geometry between them to have a well defined 
background. The logarithmic warp factor at far IR tells us that the geometry is 
influenced by one or a set of coincident D7 branes. These seven branes wrap the T^'^ 
as in branch 2 of |[13| while extending in the radial r direction and filling up four 
Minkowski directions (see eq (3.9) of ljl3[)- In particular the embedding equation for 
a D7 brane is given by |13[ 



^3/2 ^i^-^,-^,) .^ ^ Sin ^ = /i 

2 2 



(2. 



where /i is a parameter. For supersymmetric case /i would be related to the deforma- 
tion parameter of the conifold. Since we don't require supersymmetry we can take 
fi to be arbitrary^. Different values of /i will tell us how far the D7 branes are from 
the origin r = 0. The Nf D7 branes may have Nf different locations given by Nj 
different values of fi or D7 branes may be coincident with just a single value of fi. The 



^The issue of supersymmetry is a little subtle here. The susy can of course be broken by choosing 
a different fj,, but can also be broken by choosing the right fi but separating the wrapped D5 branes 
along (02, 4'2) directions. One may say that if we allow bound states of D5 and D7 branes we might 
restore zero temperature susy. Alternatively we can consider the seven-branes to be oriented as 
in which is related to our far IR configuration. In heavy fundamental quarks could still 
restore susy. In this paper we will not consider the seven-brane configurations of p8|. 
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positions of the seven branes are therefore parametrised by the coordinate z. Since 
the seven branes are in branch 2 of their positions can be precisely parametrised 
by the internal coordinates (^2? 02)- Thus the seven branes stretch along r and can 
be placed at any point on the (6'2,02) plane^. Because of this distribution, as we 
shall see shortly, axion-dilaton field runs with the coordinate z and the running is 
determined by F theory. 

Our UV cap, in the full F-theory picture, should allow a distribution of seven 
branes that could eventually reproduce the warp factors ( p.7| ). This is however not 
the only requirement: we also want to study the potential of heavy quarkonium type 
bound states in our theory (this means that we need to study the bound states of 
very heavy quark-antiquark pairs). Which in turn implies that we require a set of 
seven branes as far away from the origin as possible (or, at high temperature, as 
far away as possible from the black hole horizon). There are a few possible ways 
to distribute the seven branes that might be able to reproduce the required picture. 
The simplest way would be to distribute the seven branes as in Figure 1 below. This 



D7 Branes 




62' 4>2 



Figure 1: Simplest way to distribute localised seven branes in our model. The seven branes 
wrap the internal sphere parametrised by {9i,(j)i) and are stretched along the spacetime 
directions. Their extensions along the radial directions are parametrised by // as in embed- 
ding equation above. The coordinate rmin denote the distance of the nearest seven brane 
from the black hole horizon. A string stretched between this seven brane and the black 
hole horizon is the lightest fundamental quark in our model. The heaviest quark, on the 
other hand, will be from the seven brane that is farthest from the horizon. A string whose 
two ends lie on such a seven brane will form a quark antiquark bound state. The temporal 
evolution of such a string will determine the Wilson loop in our picture. 



^In actual case the embedding is a union of branch 1 and branch 2. Therefore the seven-branes 
will trace a complicated surface in (r, 9i, (pi, tp) plane. For simplicity we will assume the embedding 
to be given by branch 2 of Later on when we study fluxes, the non-trivial nature of the 

seven-brane embeddings will become important. 
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picture, although simple and desirable, however does not quite suffice for us because 
we need a configuration of seven branes that could interpolate between the IR and 
UV configurations. One of the simplest way to have an interpolating geometry using 



the configurations studied in |13| is to make the seven branes delocalised along the 



(r, 6'2, 02) directions and call the resulting quantity as Nf{r, 6*2, ^2)- This means that 



Nj{t) 



iV/(r,^2,02; 



sin Ori 



(2.9) 



An immediate way to realise such a configuration is given in Figure 2 below. Such a 
configuration has been advocated in some recent works (see for example [^, |2T| where 
the delocalised seven branes are embedded via the Kuperstein embeddings 1^). Such 



Distribution of seven branes 




Black Hole 



02' 4>2 



Figure 2: Complete delocalisation of the seven branes along {r,92,(p2) directions. 



a configuration of seven branes, although useful for many purposes, unfortunately 
still does not quite suffice for us because heavy quarks in such a scenario would tend 
to go to configurations of lighter quarks spontaneously. Furthermore we want to 
impose the F-theory constraint, for scales r > f: 



r>r 



24 (2.10) 



which would be a little difficult to impose in the fully delocalised scenario^. Therefore 
the configuration that we would be mostly interested in is given in Figure 3. In this 
picture, which should be viewed as a cross between the earlier two figures, every 
individual set of seven branes are delocalised a little bit. The F-theory constraint on 
the number of fiavors i.e (|2.10| ) can be easily imposed without making Nf{r,62,4>2) 



''in fact F-theory can allow number of seven-branes to be arbitrarily large. For this case we need 
to carefully study the singularity structure of the underlying manifold. Here, for most of the paper, 
we will restrict ourselves to 24 seven-branes. This means that gs could be as small as 0.042. 
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Configuration of seven branes 




02' 4>2 



Black Hole 



Figure 3: Another configuration of seven branes where the delocalisation is milder com- 
pared to the earlier picture. The local minima of every set of seven branes help us to study 
various configurations of quark antiquarks pairs. 



arbitrarily small. The final picture that we want to emphasise which would capture 
the underlying dynamics is given as Figure 4 below. The figure is a slight variant 
of the previous figure. We have divided our geometry into three regions of interest: 
Regions 1, 2 and 3. Region 1 is basically the one discussed in great details in 



Region 3 
Region 2 , '' 



Region 1 



Configuration of seven branes 



X 



.0123 





Black Hole 



02' ^1 



Figure 4: This figure, which is a slight variant of the previous figure, shows the various 
regions of interest. As should be clear, most of the seven branes lie in Region 3, except for 
a small number of coincident seven branes that dip till rmin i.e Region 1. The interpolating 
region is Region 2. The detailed backgrounds for each of these regions are given in the text. 
Note however that, although we have emphasised Region 1 more, we will only consider the 
case where Region 3 >> Region 1 + Region 2. 



In this region there is one (or a coincident set of) seven brane(s). The logarithmic 
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dependences of the warp factor and fluxes come from these coincident (or single) seven 
branes. In fact the logarithmic runnings of the gauge theory coupling constants also 
stem from these seven branes. 

Since we require UV free (or more appropriately, strongly coupled and conformal) 
the IR logarithmic runnings wouldn't be very desirable. Therefore the UV cap in the 
full F-theory framework is depicted as Region 3 in the above figure. In this region 
we expect all the seven branes to be distributed so that axio-dilaton has the right 
behavior. We also expect vanishing H^s and Hrh fields (just like the AdS cases). 

It is clear that one cannot jump from Region 1 to Region 3 abruptly. There 
should be an interpolating geometry where fluxes and the metric should have the 
necessary property of connecting the two solutions. This is Region 2 in our figure 
above. For all practical purposes, we expect Region 3 to dominate, in other words. 
Region 3 should be greater than both Regions 1 and 2 combined together. In such a 
scenario analysis of Wilson loop for heavy quark - antiquark bounds states would be 
easy: we wouldn't have to worry too much about the intermediate regions. Another 
big advantage about our UV cap is related to the issues raised in . Since the H^s 
and the axio-dilaton fields have well defined behaviors at large r, there would be no 
UV divergences of the Wilson loops in our picture! Therefore our configuration can 
not only boast of holographically renormalisability, but also of the absence of Landau 
poles and the associated UV divergences of the Wilson loops. 

In the following, let us therefore discuss the backgrounds for all the three regions 
in details. 



2.1 Region 1: Fluxes, Metric and the Coupling Constants Flow 

The background for Region 1 is discussed in details in [|13|, so we will be brief. All 
the logarithmic behaviors for the fluxes and metric come from the single set of seven 
branes. The metric has the following typical form: 



ds = —= 



gi{r)dt^ + dx^ + dy^ + dz^ 



92[r 



r^dr^ + dMl 



(2.11) 



where gi{r) are the black- hole factors and dM-i is the metric of warped resolved- 
deformed conifold (see the form in eq (3.5) of [|^). The internal space retains its 
resolved-deformed conifold form upto 0{gsNf). Beyond this order the internal space 
loses its simple form and becomes a complicated non-Kahler manifold. The warp 
factor to this order, in terms of NJ^, Meg (see eq (3.10) of |jl3[| for details), is: 



h 



1/ 

^4 



1 + 



—los,r 

2-kN ^ 



1 + 



^sNf 
27r 



logr -|- 



+ 



log 



sm — sm — 
2 2 



(2.12) 



As discussed in [jT3l , the background has all the type IIB fluxes switched on, namely, 
the three-forms, five-form and the axio- dilaton. Both Nf and M^s are different from 
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Nf and M. We will give detailed reason for this when we discuss the full geometry 
in the next two subsections. The three-form fluxes are: 



Fs = 2Meff Ai 1 



■cff 



2n ^°^") '^^2 



]- (sin 6i dOi A d(j)i - Bi sin 62 d62 A #2) 



An 



dr f 9 9 \ 

— A2 — A A ( cot — sin 6*2 d(f)2 — B2 cot — sin 9i dcpi ] 
r \ 2 2 J 



SgsM^sNf f 9o 9^ \ 

— ^Ag sin 9i sin ^2 cot — d9i + B3 cot — ^^2 A A (^02 (2.13) 

Stt \ 2 2 J 

. A , ^9sNf ^ _ QsNf , . ^1 . 92\dr 



ifs = GgsAiM^s 1 H log ^ + o sin— sin— — 

\ An 2tt 2 2 I r 



A^ I sin 9i d9i A d(j)i - B4 sin 92 d92 A 



3^?,2MeffiVf (dr ^ 

A5 ( — A - -de^ 



in 



9 9 
A I cot ^ d92 - Bg cot y d9i 



where F3 = F3 — CqH^, Cq being the ten dimensional axion and the so-called asym- 
metry factors Ai,Bi are given in eq. (3.83) of []T^ (see also JT^). The axio-dilaton 
and the five- form fluxes are: 



Co 



Nf 
An 



{ij-(f)i- 02) 

1 / 6 . 2 4N , f . ^1 . ^2 

log (r + 9a r j — ;:^log ( sm — sm 



Svr 



2n 



[d^x A dh~^ + *{d'^x A rf/i"^ 



(2.14) 



with a being the resolution parameter of the internal space that depends on the 
horizon radius r/^ as a = a{rh) + 0{g'^MeeNf). Once we consider the slice: 



9i — 92 



vr, 



0, 







(2.15) 



the background along the slice simplifies quite a bit. To 0{gsNf) the background is: 



h = 




1 + 


^9sM^s^ 
2nN,s 




= Fs = 


- Co = 






1 

9s 


Nf 
Svr 



SgsNf ( 1 
logr { 1 H I logr + - 







log (r^ + 9aV^) 



(2.16) 



alongwith F5 given by ( p.l4| ). The simplicity of the background is the reason why our 
analysis of the mass and the drag of the quark in |T3 were straightforward enough 
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to see the underlying physics, yet were not afflicted by problems like UV divergences 
of p3|^. Note that the logarithmic RG flows of the two couplings come from the 



logarithmic Bjsrs held, leading to conflnement at the far IR (at zero temperature). 
In the following, to avoid clutter, {N, Nf, M) would denote their effective values. 

2.2 Region 2: Interpolating Region and the Detailed Background 

To attach a UV cap that allows a vanishing beta function we need at least a conflgu- 
ration of vanishing NS three-form. This cannot be abruptly attached to Region 1: we 
need an interpolating region. This region, which we will call Region 2, should have 
the behavior that at the outermost boundary the three-forms vanish, while solving 
the equations of motion. The innermost boundary of Region 2 — that also forms the 
outermost boundary of Region 1 — will be determined by the scale associated with 
the mass of the lightest quark, mo, in our system. In terms of Figure 4, this is given 
by region in the local neighborhood of r^i^ = tuqTq^ + r/i, where Tq and r/i are the 
string tension and the horizon radius respectively. We have already discussed some 
aspects of this in our previous paper [|l3l when we discussed the issue of UV caps. It 



is now time to spell this in more details. 

The structure of the warp factor should be clear from [|13[. We expect the form 
to look like (|2.1| ) discussed earlier. For our purpose, it would make more sense to 
rewrite this in such a way that the radial r dependence shows up explicitly. For this 
we need to flrst deflne two functions /(r) and M(r) as (see Figure 5): 

o{r-ro) 

fir) ^ ^^^air-ro) ^ M{r) ^ M[l-f{r)], a»l (2.17) 

where the scale tq will be explained below and M is as before related to the effective 
number of flve-branes (or the RR three-form charge). Note that for r << tq, /(r) ^ 
e''~^°, whereas for r > tq, /(r) ^ 1. Thus for r smaller than the scale tq, /(r) is a 
very small quantity; whereas for r bigger than the scale tq, /(r) is identity. In terms 
of M(r) this means that for r < ro, M(r) ^ M whereas for r > ro, M(r) — )■ 0. This 
will be useful below. 

Using these functions, we see that the simplest way in which logarithmic behavior 
along the radial direction may go to inverse r behavior, is when the warp factor has 
the following form: 

h = ^i±M!l±f^y (2.18) 

a 

where Cj are constant numbers, and the denominator can be mapped to r(o) deflned 
in ( p.l|) with e(a) functions of gsNf,M,N and the resolution parameter a. L^s are 



°0n the slice (|1|) the pull-backs of the S-fields are zero. This means that Wilson loops or 
other equivalent constructions could be carried out without any interference from the logarithmic 
B-fields. 
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functions of the angular coordinates {6i, For other details see [13|. The warp 

factor h has the required logarithmic behavior as long as the exponents of r are 
small and fractional, and indeed switches to the inverse r behavior as soon as the 



exponents become integers. In [13 we gave some examples where the exponents are 
small and fractional numbers, and alluded to the case where they become integers^. 
Since Nj is a delocalised function, this behavior could be naturally realised now and 
would eventually give way to the required inverse r behavior of the warp factor in 
Region 3. Its at least clear that such a behavior of the warp factor do solve the 



background supergravity equations of motion near r = rmin (see |]13| for a concrete 
example, and we will give more details on this below), however what we want to 
know whether such a behavior of the warp factor is generically a solution to EOM, 
or we need to add sources to the theory. It will turn out that we need to add sources 
at the outermost boundary of region 2. Question now is to figure out consistently 
the specific point in the radial direction beyond which Region 3 would start. This 
way we will know exactly where to add the sources and the AdS cap. 

The demarcation point can be found easily by looking at the behavior of H^s 
and Hfiji. For this we need to use the functions ( p.l7| ) to write the RR three- form. 
Our ansatze for then is: 



3 

a„ 

2TTr9' 



^ a \ a / 



A — — > — dr A A cot — sm Oo dcho — > cot — sm 61 

a \ Q 

2^ g^^e,./ sm e, sm 62 cot -de, + }_^^ cot - dO^ 

a \ a / 



Ad0iA#2 (2.19) 



where = 1 + ^ and (c^, ha) are constants. One may also notice three things: 
first, how the internal forms get deformed near the innermost boundary of the region, 
second, how the function /(r) appears for all the components, and finally, how Nj 
is, as before, not a constant but a delocalised function^°. The function /(r) becomes 
identity for r > tq and therefore F3 — t- for r > vq. For r < tq, the corrections coming 
from /(r) is exponentially small. Integrating F3 over the topologically non-trivial 
three-cycle: 



-e^ A sin 61 dOi A - V ^ sin 62 rf^2 A (2.20) 

2 \ ^<^(a) J 

we find that the number of units of RR flux vary in the following way with respect 



^See the section on holographic renormahsability in 



^^We will soon see that Nr in fact is the effective number of seven-branes. 
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to the radial coordinate r: 



Mtot(r) = M(r) + — ^] V ^ (2.21) 

which is perfectly consistent with the RG flow, because for r < tq, and r ^ re sssM ^ 
Mtot decreases precisely as M — Nf as the correction factor e^~^° coming from /(r) is 
negligible. For r > rg, Mtot shuts off completely. This also means that below tq, the 
total colors N decrease by Mtot exactly as one would have expected for the RG flow 
with Nf flavors. Using similar deformed internal forms, one can also write down the 



For a=0.5 

For ct=0.75 

For a=l 

For ct=1.25 



Figure 5: A plot of the /(r) function for tq = 5 in appropriate units, and various choices 
of a. Observe that for large a the function quickly approaches 1 for r > vq. 



ansatze for the NS three-form. This is given as: 
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1 9a'^gsNf 



log ( sin y sin y ) + hir) 



sin 6^1 d9i A A (ir 



gs dM{r) 
127r rfr 



2 + 'dgsNf - 



1 9a^gsN f 
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with {ka,.:,qa) being constants and 6„ = Xlm T^flr where a^m = anmia^, 9sNf) 
and ?m = 'emigsNf). The way we constructed the three-forms imply that is 
closed. In fact the 0{df) terms that we added to (|2.22|) ensures that. However F3 
is not closed. We can use the non-closure of F3 to analyse sources that we need 
to add for consistency. These sources should in general be {p, q) five-branes, with 
(p, q) negative, so that they could influence both the three-forms and since the ISD 
property of the three-forms is satisfied near r = rmm the sources should be close 
to the other boundary. A simplest choice could probably just be anti five-branes 
because adding anti D5-branes would change F3, and to preserve the ISD condition, 
i^s would have to change accordingly. Furthermore, as we mentioned before, as 
r — >^ ro, both i^s = F3 — )■ 0. Therefore r = tq is where Region 2 ends and Region 3 
begins, and we can put the sources there. They could be oriented along the spacetime 
directions, located around the local neighborhood of r = tq and wrap the internal 
two-sphere (6'i,0i) so that they are parallel to the seven-branes. However, putting 
in anti D5-branes near r = Tq would imply non-trivial forces between the five-branes 
and seven-branes as well as five-branes themselves. Therefore if we keep, in general 
the (p, q) the five-branes close to say one of the seven-brane then they could get 
dissolved in the seven-brane as electric and magnetic gauge fluxes and F^^) 

respectively. Thus the seven-brane soaks in the five-brane charges, which in turn 
would mean that F3 in ( |2.19| ) and i^s in (|2.22| ) will satisfy the following EOMs: 



dF^ = A A2(2) -rf(Re r) A//3 
d*H3 = A A2(^) -t/(C4 AF3) (2.23) 

where the tension of the seven-brane is absorbed in A2{z), which is the term that 
measures the delocalisation of the seven branes (for localised seven branes this would 
be copies of the two-dimensonal delta functions) and r is the axio-dilaton that we 
will determine below. In addition to that d * F3 will satisfy its usual EOM. For all 
the analysis in this paper we will also assume: 

ko-r-minl < ai, \rmin-rh\ < a2, Region 3 >> 01 + 02 (2.24) 

to be our approximation. This way, as we said before. Region 3 will dominate our 
calculations. 

However the above set of equations ( ^^.23] ) is still not the full story. Due to the anti 



GSO projections between anti-D5 and D7-brane, there should be tachyon between 
them. It turns out that the tachyon can be removed (or made massless) by switching 
on additional electric and magnetic fluxes on D7 along, say, (r, i/j) directions! This 
would at least kill the instability due to the tachyon, although susy may not be 



restored. For details on the precise mechanism, the readers may refer to [Q. But 
switching on gauge fluxes on D7 would generate extra D5 charges and switching on 
gauge fluxes on anti-D5s will generate extra D3 charges. This is one reason why we 
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write (A^, Nf, M) as effective charges. This way a stable system of anti-D5s and D7 
could be constructed. 

To complete the rest of the story we need the axio-dilaton r and the five-form. 
The five-form is easy to determine from the warp factor h ( glsD using (|2.14|) . The 



total five-form charge should have contribution from the gauge fluxes also, which in 
turn would effect the warp factor. For regions close to rmm it is clear that r goes as 
^-gsNf -^Y^QYe z is the embedding (|2.8| ). More generically and for the whole of Region 
2, looking at the warp factor and the three-form fluxes, we expect the axio-dilaton 



to go as^"*^: 

r = [bo + hf{r)]J2^ (2.25) 

a 

where bi are constants and Ca are functions of the internal coordinates and are 
complex. These Ca and the constants 6j are determined from the dilaton equation 
of motion fEl, |25ll : 



.2 _ Vr ■ Vr 4fi;2Q(Im 63^7 



r = — ' + (p, q) sources (2.26) 

zim r y~g or 

where tilde denote the unwarped internal metric gmn, and is the action for the 
delocalised seven branes. The /(r) term in the axio-dilaton come from the {p,q) 
sources that are absorbed as gauge fluxes on the seven-branes^^ . Because of this 
behavior of axio-dilaton we don't expect the unwarped metric to remain Ricci-flat 
to the lowest order in QsNf. The Ricci tensor becomes: 

~ _ 2 9(mdn)T 2 f fD7 _ 1~ fDj] , ^2 ( fip,q}5-hrs.ne _ 1~ f;(p,g)5-brane 



(2.27) 



where we see that T^^r picks up terms proportional to e^^.^ and derivatives of /(r), Nf{r) 
implying that to zeroth order in QsNf the interpolating region may not remain Ricci- 
flat. However since the coefficients are small, the deviation from Ricci-flatness is 
consequently small. In this paper we will not give the explicit form for Ca, La etc 
but it should be clear from our above discussions that EOMs are easily satisfied. The 



11 



One may use this value of axio-dilaton and the three-form NS fluxes (2.22) to determine the 



beta function from the relations (2.2). To lowest order in gsNf we will reproduce the SV beta 
function (2.3) as expected. Notice that for r > ro the beta function does not vanish and both 
the gauge groups flow at the same rate. This will be crucial for our discussion in the following 
subsection. 

^^The r"*^!") behavior stems from additional anti seven-branes that we need to add to the existing 
system to allow for the required UV behavior from the F-theory completion. The full picture will 
become clearer in the next sub-section when we analyse the system in Region 3. 
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one last thing to check would be the equation for the warp factor. This is given by 
the five-form equation of motion: 



d * dh^^ = if3 A F3 + tr (f(^) A F^^^ - 7^ A 7^) A2(z) + k^q tr F^'^^A^iS) 

(2.28) 

where F^^^ is the seven-brane gauge fields that we discussed earlier, F^^^ is the (p, q) 
five-brane gauge fields required for the proper interpretation of the colors in the gauge 
theory side^^, TZ is the pull-back of the Riemann two- form, and A^IS) is the term 
that measures the delocalisation of the dissolved {p, q) five-branes over the space S 
embedded in the seven-brane (again for localised five-branes there would be copies 



of four- dimensional delta functions). The H^, A F3 term in ( |2.28|) is proportional to 

MHr) 



This is precisely the form for the warp factor ansatze (|2.18| ) with the / (r) 
term there accounting for the M^(r) term above. This way with the warp factor 
( p.l8| ) and the three-forms ( p.l9| ) and ( p.22| ) we can satisfy ( p.28| ) by switching on 



small gauge fluxes on the seven-branes and five-branes. 

Therefore combining ( |2.18| ), ( |2.19| ), ( |2.22| ), ( |2.25| ) and the five- form, we can pretty 
much determine the supergravity background for the interpolating region r-ann <r < 
tq. At the outermost boundary of Region 2 we therefore only have the metric and 
the axio-dilaton. Both the three-forms exponentially decay away fast, giving us a 
way to attach an AdS cap there. 

2.3 Region 3: Seven Branes, F-Theory and UV Completions 

The interpolating region. Region 2, that we derived above can be interpreted alter- 
natively as the deformation of the neighboring geometry once we attach an AdS cap 
to the OKS-BH geometry. The OKS-BH geometry is the range < r < rmin and 
the AdS cap is the range r > tq. The geometry in the range rmin < t < tq is the 
deformation. Such deformations should be expected for all other UV caps advocated 
in [|l3l. In this section we will complete the rest of the picture by elucidating the 
background from r > tq in the AdS cap. But before that let us give a brief gauge 
theory interpretation of background'^^. 

For the UV region r > tq we expect the dual gauge theory to be SU {N + M) x 
SU{N + M) with fundamental flavors coming from the seven-branes. This is because 
addition of (p, q) branes at the junction, or more appropriately anti five-branes at 
the junction with gauge fluxes on its world-volume, tell us that the number of three- 
branes degrees of freedom are + M, with the M factor coming from five-branes 



-'^^In fact one should view the gauge fluxes on the seven-branes and the five-branes as the total 
gauge fluxes that are needed to stabilise the system. We will see in the next subsection that the 
full stabisation would require additional fluxes, but the structure would remain the same. 

^*The discussion in the following paragraph is motivated by a correspondence that we had with 
Peter Ouyang. We thank him for his comments. 
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anti-five-branes pairs. Furthermore, the SU{N + M) x SU{N + M) gauge theory 
will tell us that the gravity dual is approximately AdS, but has RG flows because 
of the fundamental flavors (This RG flow is the remnant of the flow that we saw in 
the previous subsection. We will determine this in more details below). At the scale 
r = ro we expect one of the gauge group to be Higgsed, so that we are left with 
SU{N + M) X SU{N). Now both the gauge fields flow at different rates and give 
rise to the cascade that is slowed down by the Nf flavors. In the end, at far IR, we 
expect confinement at zero temperature. 

The few tests that we did above, namely, (a) the flow of and M colors, (b) 
the RG flows, (c) the decay of the three-forms, and (d) the behavior of the dual 
gravity background, all point to the gauge theory interpretation that we gave above. 
What we haven't been able to demonstrate is the precise Higgsing that takes us to 
the cascading picture. From the gravity side its clear how this could be interpreted. 
From the gauge theory side it would be interesting to demonstrate this. 

Coming back to the analysis of Region 3, we see that in the region r > rg we do 
not expect three-forms but we do expect non-zero axio-dilaton. These non-zero axio- 



dilaton come from the rest of the seven branes. As mentioned in the complete 



set of seven-branes should be determined from the F-theory picture [27| to capture 



the full non-perturbative corrections. This is now subtle because the seven-branes 



are embedded non-trivially here (see (|2.8|) ). A two-dimensional base, parametrised 



by a complex coordinate z, on which we can have a torus fibration: 

= + xF{z) + G{z) (2.29) 



can be identified with the z coordinate of ( |2.^ ). This way vanishing discriminant 
A of ( p.29| ) i.e A = AF^ + 27G^ = 0, will specify the positions of the seven-branes 



exactly as ( |2.8|) . Here we have taken F{ degree eight polynomial in z and G{z) 

as a degree 12 polynomial in z. The delocalisation Nf{r,62,(f)2) should be thought 
of somewhat as the distribution of bunches of seven branes along (^2,02) directions 
with varying sizes along the radial r direction such that ( |2.10| ) is maintained with 



the deviation 5 = f — ro a finite but not very large number. 

As is well known, embedding of seven-branes in F-theory also tells us that we can 
have SL{2, Z) jumps of the axio-dilaton. We can define the axio-dilaton r = Co+ie~'^ 
as the modular parameter of a torus fibered over the base parametrised by the 
coordinate z. The holomorphic map^^ from the fundamental domain of the torus to 
the complex plane is given by the famous j-function: 

. ^ [9g(r) + e^(r) + e|(r)f ^ 4(24F(z))3 



15 



Holomorphic in r, the modular parameter. 
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where Qi,i = 1,2,3 are the well known Jacobi Theta-f unctions and rj is the Dedekind 
r^-f unction: 



r]{T) = giJJ(l-g"), q 



(2.31) 



For our purpose, we can write the discriminant A (2;) and the polynomial F{^z) gener- 
ically as: 



A(z) = 4F3 + 27G^ 



24 

n 



'3)1 



F{z) = fefjl 



Z - Zi 



(2.32) 



so that when we have weak type IIB coupling i.e t = Cq + ioo, ji^r) ^ e and 



using ( p.30|) the modular parameter can be mapped to the embedding coordinate z 
as: 



+ ^ log (55926ar^) - 



n=0 



n=l 



\i=i i=i 



(2.33) 



where C„ = Cn{Oi, (l)i,ip) and = 'Dn{9i,(l)i,'ip) are real functions and r = r^/^. To 
avoid cluttering of formulae, we will use r instead of r henceforth unless mentioned 
otherwise. So the coordinate r will parametrise Region 3, and t = XI '^"t^" ■ 

The above computation was done assuming that z > {zi,Zj), which at this stage 
can be guaranteed if we take 6'i 2 small. This gives rise to special set of configurations 
of seven-branes where they are distributed along other angular directions. However 
one might get a little worried if there exists some zj = Zo related to the farthest 
seven-brane(s) where the above approximation fails to hold. This can potentially 
happen when we try to compute the mass of the heaviest quark in our theory. The 
question is whether we can still use the r derived in (|2.33|) , or we need to modify the 
whole picture. 

Before we go into answering this question, the choice of z bigger than {zi,Zj) 
already needs more convincing elaboration because allowing 6'i 2 small is a rather 
naive argument. The situation at hand is more subtle than that and, as we will 
argue below, the picture that we have right now is incomplete. 

To get the full picture, observe first that z being given by our embedding equation 
( p.8|) , means that if we want to be in Region 3, we need to specify the condition 
r > To in the defination of z. This way a given z will always imply points in Region 3 
for varying choices of the angular coordinates {di,(j)i,ip). However similar argument 
cannot be given for any choices of {zi,Zj). A particular choice of {zi,Zj) may imply 
very large r with small angular choices or small r with large angular choices. Thus 
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analysing the system only in terms of the r coordinate is tricky. In terms of the full 
complex coordinates, z > {zi,Jj) would mean that we are always looking at points 
away from the surfaces given hj z = Zi and z = 'zj. 

What happens when we touch the z = Zi surfaces? For these cases F{zi) — ?■ 
and therefore we are no longer in the weak coupling regime. For all F{zi) = imply 
j(r) —7- which in turn means r = exp («7r/3) on these surfaces. These are the 



constant coupling regimes of |28| where the string couplings on these surfaces are 
not weak. On the other hand, near any one of the seven-branes z = Zj we are in the 
weak coupling regimes and so (|2.33|) will imply 



t{z) = — log {z - Zj) ioo (2.34) 

which of course is expected but nevertheless problematic for us. This is because we 
need logarithmic behavior of axio-dilaton in Region 2, but not in Region 3. For a 



good UV behavior, we need axio-dilaton to behave like ( 2.33 ) everywhere in Region 
3. 

In addition to that there is also the issue of the heaviest quarks creating additonal 
log divergences that we mentioned earlier. These seven branes are located aX z = 
Zj = Zo, and therefore if we can make the axio-dilaton independent of the coordinates 
Zo then at least we won't get any divergences from these seven-branes. It turns out 
that there are configurations (or rearrangements) of seven-brane(s) that allow us to 
do exactly that. To see one such configuration, let us define F{z),G{z) and A(z) in 
( p.32|) in the following way: 

7 10 

F{z) = {z-Zo)l[iz-Zi), G{z) = {z-Zoyl[{z-z,) 

1=1 1=1 

21 

A(^) = {z~Zo)''l[{z-I,) (2.35) 

i=i 

which means that we are stacking a bunch of three seven-branes at the point z = Zo, 
and 

21 7 10 

= 4Y[{z-z,f + 27{z-Zo)\[{z-hf (2.36) 

j=l i=l i=l 

implying that the axio-dilaton r becomes independent of io and behaves exactly as 
in ( p.33| ) with (z, j) in ( |2.33| ) varying upto (7, 21) respectively. 



The situation is now getting better. We have managed to control a subset of 
log divergences. To get rid of the other set of log divergences that appear on the 
remaining twenty-one surfaces, one possible way would be to modify the embedding 
). Recall that our configuration is non-supersymmetric and therefore we are not 



- 21 - 



required to use the embedding (2^). In fact a change in the embedding equation will 



also explain the axio-dilaton choice (|2.25|) of Region 2. To change the embedding 



equation ( |2.8| ) we will use similar trick that we used to kill off the three-form fluxes, 
namely, attach anti-branes. These anti seven-branes^^ are embedded via the following 
equation: 

^3/2giW-<Ai-02)sin I Sin I = roe^® (2.37) 

where is some angular parameter, and could vary for different anti seven-branes. 
The above embedding will imply that their overlaps with the corresponding seven- 
branes are only partial^^. And since we require 

it will appear effectively that we can only have seven-branes in Regions 1 and 2, 
and bound states of seven-branes and anti seven-branes in Region 3.^*^ This way the 
axio-dilaton in Region 3 will indeed behave as ( |2.33| ) for all z (except for the above 
mentioned seven points). 

There are two loose ends that we need to tie up to complete this side of the 
story. The first one is the issue of Gauss' law, or more appropriately, charge con- 
servation. The original configuration of 24 seven branes had zero global charge, but 
now with the addition of anti seven-branes charge conservation seems to be prob- 
lematic. There are a few ways to resolve this issue. First, we can asume that that 



branes wrap topologically trivial cycles, much like the ones of Then charge 

conservation is automatic. The second alternative is to isolate six seven-branes using 
some appropriate F and G functions, so that they are charge neutral. This is of 



course one part of the constant coupling scenario of p9|. Now if we make the {62, 02 



directions non-compact then we can put in a configuration of 18 seven-branes and 
anti seven-branes pairs together using the embeddings ( p.8| ) and ( |2.37| ) respectively. 
The system would look effectively like what we discussed above. Since the whole 
system is now charge neutral, compactification shouldn't be an issue here. 

The second loose end is the issue of tachyons between the seven-brane and anti 
seven-brane pairs. Again, as for the anti-D5 branes and D7-brane case [^, switching 
on appropriate electric and magnetic fluxes will make the tachyon massless! Therefore 
the system will be stable and would behave exactly as we wanted, namely, the axio- 
dilaton will not have the log divergences over any slices in Region 3. 



-'^^They involve both local and non-local anti seven-branes. 

"'^''For example if we have a seven-brane at z = zi such that lowest point of the seven brane is 
r = |zip/'^ < To, then the corresponding anti-brane has only partial overlap with this. 

"'^^Of course this effective descrition is only in terms of the axio-dilaton charges. In terms of the 
embedding equation for the seven-branes ( |2.8| ) this would imply that we can define z with r > 
and 2j with r < r^. 
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This behavior of axio-dilaton justifies the r~^('^) in (|2.25|) in Region 2. So the full 
picture would be a set of seven-branes with electric and magnetic fluxes embedded 
via ( p.8|) and another set of anti seven-branes embedded via ( p.37|) lying completely 
in Region 3. 

Thus in Region 3 both the three-forms vanish and therefore gi = g2 = gyu with 
gi,g2 being the couplings for SU{N + M), SU {N + M). From ( p^ ) we can compute 
the /3-function for gYM as: 

^^^^^^ = 91^ - T6^ ^ (2-38) 

n=l 

where A is the usual RG scale related to the radial coordinate in the supergravity 
approximation. For A — oo, /3((7ym) implying a conformal theory in the far 
UV. We can fix the 't Hooft coupling to be strong to allow for the supergravity 
approximation to hold consistently at least for all points away from the z = Zi,i = 
1, 7 surfaces. 

Existence of axio-dilaton r of the form ( p.33| ) and the seven-brane sources will 



tell us, from (|2.27|) , that the unwarped metric may not remain Ricci flat. For example 



it is easy to see that 



. ^ f + ^ o (L) (2.39) 

n,m=l ^ ^ 

where the last term should come from the seven-brane sources and, because of these 
sources, we don't expect Tirr to vanish to lowest order in g^Nj}^ The term At> is 
given by the following infinite series: 



k,l=l k,l,p,q=l 

Similarly one can show that 

/^afe - ^ 2^ + ^ I J ^ ' 

for (a, b) 7^ r. For TZrb similar inverse r dependence can be worked out. In the far 
UV we expect the unwarped curvatures should be equal to the AdS curvatures. The 
warp factor h on the other hand can be determined from the following variant of 
O: 



d * dh-^ = tr (F^^) A F^^) - 7^ A 7^) A2iz) + ... (2.42) 
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Although, as discussed before, the deviatfon from Ricci flatness will be very small. 
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because we expect no non-zero three- forms in Region 3. The dotted terms are the 
non-abehan corrections from the seven-branes. As r is increased i.e r » tq, we 
expect F^^^ to fall-off (recall that they appear from the anti (1,1) five-branes located 
in the neighborhood of r = tq) and therefore can be absorbed in TZ. Once we embed 
the seven-brane gauge connection in some part of spin-connection, we expect 

nh-^ = o 

from the non-abelian corrections via pull-backs. Solving this will reproduce the 
generic form for h: 

with a constant and a^'s are suppressed by powers of QgNf. More details on this 
is given in the Appendix A and B. At far UV we recover the AdS picture implying 
a strongly coupled conformal behavior in the dual gauge theory. 

From the above discussion we can conclude that the warp factor and the axio- 
dilaton will have the inverse r behavior. We will use this background to do the 
Wilson loop computation in the next section. 




ai{ip,9i,(j)i) 



i=l 



(2.44) 



3. Heavy Quark Potential from Gravity 

Before we go into the actual computation of the Wilson loop, let us point out some 
generic standard arguments that map the Wilson loop computation to the string 
action and then to the quark anti-quark potential. 

Consider the Wilson loop of a rectangular path C with spacelike width d and 
timelike length T. The timelike paths can be thought of as world lines of pair of 
quarks QQ separated by a spatial distance d. Studying the expectation value of the 
Wilson loop in the limit T — )■ oo, one can show that it behaves as 

(WiC)) ~ expi-TEQQ) (3.1) 

where Eqq is the energy of the QQ pair which we can identify with their potential 
energy Vqq(c?) as the quarks are static. At this point we can use the principle of 
holography P |]5D[ and identify the expectation value of the Wilson loop with 



the exponential of the renormalised Nambu-Goto action, 

{WiC)) ~ exp(-^-a) (3.2) 

with the understanding that C is now the boundary of string world sheet. Note that 
we are computing Wilson loop of gauge theory living on flat four dimensional space- 
time x*^'^'^'^. Whereas the string worldsheet is embedded in curved five-dimensional 
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manifold with coordinates x'^'^'^'^ and r. We will identify the five- dimensional mani- 
fold with Region 3 that we discussed above. 

To be consistent with the recipe in |T0[, we need to make sure that the induced 



four dimensional metric at the boundary of the string world sheet C is flat. For an 
AdS space, this is guranteed as long as the world sheet ends on boundary of AdS 
space where the induced four dimensional metric can indeed be written as r]^^. Using 
the geometry constructed in the previous section for Region 3, we see that the metric 
is asymptotically AdS and therefore induces a fiat Minkowski metric at the boundary 
via: 

lim u'^g^^ = 7]^^ (3.3) 

n->-0 

where u = and g^i, is the full metric (including the warp factor) in Region 3. 
Thus we can make the identification ( p.2|) . Once this subtlety is resolved, comparing 
( p.l|) and ( p.2|) we can read off the potential 

cren 

VOO = # (3.4) 

Thus knowing the renormalised string world sheet action, we can compute Vqq for 
a strongly coupled gauge theory. 

The above discussion was all for gauge theory at zero temperature. What hap- 
pens when we allow non-zero temperatures? Does the above identification (|3.4| ) 
between the quark anti-quark potential and the renormalised Nambu-Goto action go 
through again? 

The answer is yes, but the derivation is a little more subtle than what we pre- 
sented for the zero temperature case. At high temperatures and density we expect 
the medium effects to screen the interaction between the heavy quark and anti-quark 
pairs. The resulting effective potential between the quark anti-quark pairs separated 
by a distance d at temperature T can then be expressed succinctly in terms of the 
free energy F{d,T), which generically takes the following form: 

Fid, T) = ad Md, T) - ^Ud, T) (3.5) 

where a is the string tension, a is the gauge coupling and fc and fs are the screening 
functions^° (see for example |^ and references therein). For the quark and the anti- 



quark pair kept at +| and — | we expect the Wilson lines W (±f) to be related to 
the free energy via: 



exp 



F{d, T) 
T 



{W^ {+i)W {-^)) ^^^^^ 



{+i)){w{-i)) 



^°We expect the screening functions fs,fc to equal identity when the temperature goes to zero. 
This gives the zero temperature CorneU potentiaL 
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In terms of Wilson loop, the free energy ( |3.5|) is now related to the renormalised 
Nambu-Goto action for the string on a background with a black-hole^^ . One may 
also note that the theory we get is a four-dimensional theory compactified on a circle 
in Euclideanised version and not a three-dimensional theory. 

3.1 Computing the Nambu-Goto Action: Zero Temperature 

Our first attempt to compute the NG action would be to consider the zero tempera- 



ture case. This means that we make the black-hole factors gi in (|2.11|) to be identity. 
The string configuration that we will take to do the required computation is given 
below in Figure 6. Note that we have configured our geometry such that the string 
is exclusively in Region 3. We will provide a stronger motivation for this soon. For 
the time being observe that the configuration in Figure 6 has one distinct advantage 
over all other configurations studied in the literature, namely, that because of the 
absence of three-forms in Region 3 we will not have the UV divergence of the Wilson 
loop attributed to the logarithmically varying B field [^. In fact even if the string 
enters Regions 2 and 1 we will not encounter any problems because there are no UV 
three-forms in our model. 

Since the system is not dynamical, the world line for the static QQ can be chosen to 
be 

= x^ = x^ = (3.7) 

and using u = 1 /r we can rewrite the metric in Region 3 as^^ : 

ds^ = g^udX^'dX" = An{4j,ei,(f)i)u''-^[-g{u)dt^ + d'^^] 



where An are the coefficients that can be extracted from the Oj in (|2.44|) , the black 
hole factor g{u) = 1 for the zero-temperature case, and ds'j^^ is the metric of the 
internal space that includes the corrections given in (|2.41|) . This can be made precise 
as 



AnU 



aiu ( 3a? 02 



2 



(3.9) 



giving us ^0 = ^^^i = ~Jl^j^2 = — and so on. Note that since Oj, 

i > 1 are of 0{gsNf) and oc ^ gsN , all Ai are very small. The r~" corrections 
along the radial direction given in ( 2.39| ) are accomodated above via BiU^ series. 



^^There is a big literature on the subject where quark anti-quark potential has been computed 
using various different approaches like pNRQCD |32), hard wall AdS/CFT [|3[ |34) and other tech- 



niques [gO|, Its reassuring to note that the results that we get using our newly constructed 
background matches very well with the results presented in the above references. This tells us that 
despite the large N nature there is an underlying universal behavior of the confining potential. 
^^We will be using the Einstein summation convention henceforth unless mentioned otherwise. 
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A B 



r =co 



C D 




A 



-0123 



^min r - a 



Figure 6: The string configuration that we will use to evaluate the Wilson loop in the 
dual gauge theory. The line A determines the actual boundary, with the line B denoting 
the extent of the seven brane. We will assume that line B is very close to the line A. The 
line C at r = To denotes the boundary between Region 3 and Region 2. Region 2 is the 
interpolating region that ends at r = r^i^. At the far IR the geometry is cut-off at r = a 
from the blown- up 5'^. As discussed in the text, the string has a maximum dip that will 
eventually lead to the confining potential between the heavy quark and the antiquark. 



Now suppose : (a, r) — )■ {x^^^^,u,'4',(j)i,9i) is a mapping from string world 
sheet to space-time. Choosing a parametization t = = t,a = = x with the 
boundary of the world sheet embedding being the path C, we see that we can have 

X° = t, = X, X"^ = X^ = 0, = X^ = i> = Q 

(X^X5) = (^1,01) = (7r/2,0), = (^2,02) = (7r/2,0) (3.10) 

which is almost like the slice ( p.l5|) that we choose in The advantage of such a 
choice is to get rid of the ackward angular variables that appear for our background 
so that we will have only a r (or u) dependent background like in (|2.16|) discussed 
before. We will also impose the boundary condition 

u{±d/2) = ^ (3.11) 

where denote the position of the seven brane closest to the boundary. The Nambu- 
Goto action for the string connecting this seven brane is: 

^•string = ^ / dadT\yl-dei [{g^, + d^cpd,cP)daX^^d,X'^] + ^£'^"5,, + J(0) 

+ e'^'daX^d.x'^ e r^r'^'"=-T„ e Fa,,.... + o(e^)] (3.12) 

where a,b = 1,2, di = 82 = The other fields appearing in the action are 
the pull backs of the NS B field Bat, the dilaton coupling J(0) and the RR field 
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strengths Fabc- Its clear that once we switch off the fermions i.e = = the 
RR fields decouple. The B^s held do couple to the fundamental string but as we 
discussed before, in Region 3 we don't expect to see any three-form field strengths. 
This is because the amount of 5ns that could leak out from Region 2 to Region 3 is: 



5ns = MS[1 - /(r)] = MS e-^^^-^'o), r > tq 
where S is the two-form: 



(3.13) 



S = gs (bi (r ) cot ^ dOi + 62 (r) cot ^ 0162] A - ^4^!^ sm ^2 ^^2 A #2 (3. 14) 
\ / 11 IZn 



+ 



3gs 

An 
9s 



1 + 9sNf - + ^2 ) log ( y sm - ) + 63 (r) 



sin 6^1 d6'i A 



127r 



1 



j.l6gaNf 



+ 1 sm ^2 c^^2 A ucp2 



^2gsNf 



and 6„ have been defined before. We see that not only 5ns has an inverse r fall off, 
but also has a strong exponential decay because a >> 1. This is the main reason 
why there are no NS or RR three-forms in Region 3, making our computation of the 
Wilson loop relatively easier compared to the pure Klebanov-Strassler case. 

On the other hand dilaton will couple additionally via the J((/>) term. Although 
this coupling of (j) is not to the X^, we can still control this coupling by arranging 
the other seven-branes such that: 



> i=l i=l 



< e 



for < n < mo 



(3.15) 



with e very small and a sufficiently big number. Under this condition the dilaton 
will be essentially constant and the axio-dilaton r would behave as: 



(3.16) 



n=l 



n>r?io 



SO that its contribution to NG action can be ignored although the Biu'' contribution 
still dominates, because the seven-branes continue to affect the geometry from their 
energy-momentum tensors and the axion charges. In this limit both string and 
Einstein frame metrics are identical and the background dilaton is 



log gs - gsVn^ 



-.U 



n+nio 



+ 0{gl) 



(3.17) 



which, in the limit (7^ — ?■ 0, will be dominated by the constant term (note that mo is 
fixed). Because of this form, the NG string will see a slightly different background 
metric as evident from ( |3.12|) . 
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Thus once the dust settles, using the metric ( p.8|) with the embedding given 
by (|3.10|) , one can easily show that at zero temperature the NG action is given by: 



5, 



NG 



T 
2^ 



^2 dx 

d 
2 



dx J 
(3.18) 



where we have used J dt = T/Tq = T (with Tq = 1 henceforth), Vn+mo = {n + 
mo)T>n+rno] ^"^^ An, Bn and Vn+mo ^rc uow defined for choices of the angular coordi- 
nates given in ( p.lOp . The above action can be condensed if we redefine: 

B„,u"' + 2g'^Vn+rn.Vi+n.Aku''''''-'-''""' + Oigt) = Giu' (3.19) 

which would mean that the constraint equation i.e c^iTj'^ = 0, being the stress- 
tensor, for u{x) derived from the action (|3.18|) using (|3.19D can be written as 




{An U^) 







implying that: 



(A^ U"^y + Gn^U"^ {^) 



u^sJ{AmU^f + g^M" u'{xy 



(3.20) 



Co 



(3.21) 



where Co is a constant, and u'{x) = This constant Co can be determined in the 
following way: as we have the endpoints of the string at x = ±d/2, by symmetry 
the string will be U shaped and if Wmax is the maximum value of m, we can define 
m(0) = Mmax and u'{x = 0) = 0. Plugging this in ( p.21| ) we get: 



Co 



•^n "-max 



(3.22) 



Once we have Co, we can use (|3.20|) to get the following simple differential equation: 



du 
dx 



± 



1 



(AnU^) 



CoyJQniU^' 

which in turn can be used to write x{u) as: 



- Cl {Ar^v!' 



1/2 



(3.23) 



x{u) = Co 



where we have used x(ur, 



dw- 



{AnW' 



1 - 



[An.W' 



-1/2 



(3.24) 



0. Now using the boundary condition given in ( ^.11| ) 



i.e x[u = u. 



7/ 



2Ur 



d/2, and defining w = Uma.xV, 



we have 



dv 



( A 7/"* if' 



A ?/" 

•^n "-max 
•^m"'max'^ 



-1/2 



(3.25) 
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At this stage we can assume all An > 0. This is because for An > we can 
clearly have degrees of freedom in the gauge theory growing towards UV, which is 
an expected property of models with RG flows. Of course this is done to simplify 
the subsequent analyses. Keeping An arbitrary will also allow us to derive the linear 
confinement behavior, but this case will require a more careful analysis. We will leave 
this for future works. Note also that similar behavior is seen for the the Klebanov- 
Strassler model, and we have already discussed how degrees of freedom run in Regions 
2 and 3. Another obvious condition is that d, which is the distance between the 
quarks, cannot be imaginary. From ( p.25| ) we can see that the integral becomes 
complex for 



'^ri "max 
•^m"max^ 



> 1 



(3.26) 



whereas for J-'{v) 
must have 



1 the integral becomes singular. Then for d to be always real we 



J^iv) < 1 



(3.27) 

We can now use, without loss of generality, = 1 and = in units of L^. Such 
a choice is of course consistent with supergravity solution for our background (as 
evident from (|3.9|) ). Therefore analyzing the condition (|3.27|) , one easily finds that 
we must have 



^(m + l)Am+3 u 



m+3 
max 



< 1 



(3.28) 



for d to be real. This condition puts an upper bound on Mmax and we can use this to 
constrain the fundamental string to lie completely in Region 3 as depicted in Figure 
5 earlier. Observe that for AdS spaces. An = for n > and hence there is no upper 
bound for Mmax- This is also the main reason why we see confinement using our 
background but not from the AdS backgrounds. Furthermore one might mistakenly 
think that generic Klebanov-Strassler background should show confinement because 
the space is physically cut-off due to the presence of a blown- up S^. Although such 
a scenario implies a Wmax for the fundamental string, this doesn't naturally lead to 
confinement because due to the presence of logarithmically varying Sns fields there 
are UV divergences of the Wilson loop. These divergences cannot be removed by 
simple regularization schemes p3| . 

Coming back to (|3.18|) we see that it can be further simplified. Using (|3.21| ), 
( |3.22|) and (|3.23| ) in (|3.18|) , we can write it as an integral over u: 

-1/2 

I J _ 11, 

1 - 



NG 



T 

IT 



du 



^AniU 



T 



iiax J e„ 



TC U 



A ?/" 

•^n "max 
•^m"'max'^ 



-1/2 



(3.29) 
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where in the second equahty we have taken v = u/umax- 

This simphfied action ( p.29D is not the full story. It is also divergent in the limit 
Co — > 0. We isolate the divergent part of the above integral ( |3.29| ) by first computing 



it as a function of €(,. The result is 



'S'ng — + S^G — 
T 1 dv 



T 1 



71 U 



max J to 



+ 



vr u 



max J eo 



A ?y" 

""mi 



-1/2 



(3.30) 



Now by expanding a/ ^mMSax^™ 



' we can compute the first integral to be 



NG 



T 1 



71 Ur 



(3.31) 



where Qo = ^o, Qi = l^i^max and so on. The second integral becomes 



oil 



T 1 



I —\fQ'r, 

Jo ^ 



1 -t; 



4 / -^""^max 

A II™- if' 



:1 -1/2 



(3.32) 



where the eo dependence here appears to O(e^); and we have set = without 
loss of generality. Now combining the result in (|3.31|) and ( p.32|) , we can obtain the 
renormalized action by subtracting the divergent term 0{l/t) in the limit eo — !■ 
and obtain the following result 



Qren 
^NG 



T 1 



1=2 



+ 



m"'max'^ 



-^w "max 

m n.m 

rrt "-max ^ 



A u 



^v™ + 0{gl) 

-1/2 



0{eo) ) (3.33) 



where the third term in ( |3.33| ), including the 0{g1) correction, is related to the action 
for a straight string in this background in the limit Qs 0. Our subtraction scheme 
is more involved because the straight string sees a complicated metric due to the 
background dilaton and non-Ricci flat unwarped metric. This effect is independent 
of any choice of the warp factor. We expect this action to be finite in the limit 

Once we have the action, we should use this to compute the QQ potential through 
( p.4| ) . Looking at ( p.25[ ) we observe that the relation between d and Wmax is parametric 
and can be quite involved depending on the coefficients An- If we have An = 0, ^„ = 
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for > 0, we recover the well known AdS result, namely: d ~ Wmax and Vqq ~ ^. 
But in general (p.25|) and ( p. 331) should be solved together to obtain the potential. 

As it stands, (pl) and (1031) are both rather involved. So to find some correla- 
tion between them we need to go to the limiting behavior of Wmax- Therefore in the 
following, we will study the behaviour of d and S^q for the cases where Wmax is large 
and small. 



3.1.1 Quark- Antiquark potential for small Wmax 

Let us first consider the case where Umax is small. In this limit we can impose 
= eUmax, SO that eo = e in all the above integrals and consequently their lower 
limits will be independent of Wmax- We can also approximate 



A^^ma^ = A + A^max = 1 + ^ 



(3.34) 



where Aq = 1 and ^2Wmax — V- Using this we can write both ( p.25| ) and ( p.33|) as 
Taylor series in 77 around 77 = 0. The result is 



d = a/t/ [do + 0.1?] + C(?7^ 
60 + biT] + 0{ri'^ 



cren 



T 



TT 



(3.35) 



with oo, cti, bo, bi are defined in the following way: 



ao 
ai 

bo 

bi 



/A 
2 



dv- 



1.1981 



2 JO 



dv- 



2 JO VT^' 



VA2 



+ 



2A2 



'A2 

1 



-1 + / dv 



1 - 



V 



-0.62 a/ A 



2aM2 



^ \2A2V^ + Q2v\\ + i;2)(i - v/r^) 







t;2(l +t.2)yr^ 



(3.36) 



where we have taken = \ and ^1 = without loss of generality. Note that all 
the above integrals are independent of 77 (or Wmax) because all 0{t) corrections are 
independent of Umax- Note also that 60 = — |^o|- In this limit clearly increasing 77 
increases rf, the distance between the quarks. For small rj, d = ao^, and therefore 
the Nambu-Goto action will become: 



cren 
'-'NG 



T 



ao\bo 



TT 



1 
d 



d + C(rf^ 



(3.37) 
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where all the constants have been defined in (|3.36| ). Using (|3.4|) we can determine 
the short- distance potential to be (recall Tq = 1): 



TT 

0.236 



'3\ 



d 



+ (0.174^2 + 0.095^2) + Oid-^) (3.38) 



which is dominated by the inverse d behavior, i.e the expected Coulombic behavior. 
Note that the coefficient of the Coulomb term is independent of the warp factor and 
therefore should be universal. This result, in appropriate units, is of the same order 
of magnitude as the real Coulombic term for the Charmonium spectra pT|, BR p2], R| 



This prediction, along with the overall minus sign, should be regarded as a success of 
our model (see also where somewhat similar results have been derived in a string 
theory inspired model). The second term on the other hand is model dependent, and 
vanishes in the pure AdS background. 

Note also that the above computations are valid for infinitely massive quark- 
antiquark pair. For lighter quarks, we expect the results to differ. It would be 
interesting to compare these results with the ones where quarks are much lighter. 

3.1.2 Quark- Antiquark potential for large Wmax 

To analyse the quark-antiquark potential for large Wmax we ffrst define a quantity 
called 2;max = ^"^^ which would be our small tunable parameter. We note that the 
smallest z^ax will come from the following equality: 



m+l)^ = 1 (3.39) 



which is the upper bound on the inequality ( |3.28|) . Furthermore since we demanded 
An > 0, the above condition will imply that the coefficients An has to quickly become 
very small because in the hmit Zmax < 1 

TTl ~t~ 1 

hm ^ 00 (3.40) 

max 



which is perfectly consistent with ( p. 9]) because higher An are proportional to higher 
powers of QsNf and therefore strongly suppressed in the limit QgNf — 0. This will 
also mean that we can retain only few of the AnS to study the potential for small 
Zjnax- In fact we will soon give an estimate of the largest n that we should keep in 
our analysis. 

To determine the distance between the quark we can again use (|3.24|) . However 



we have to be careful about a few subtleties that appear due to our choice of the scale 
^max- First of all note that we will use w = z^^g^^v in (|3.24| ). This will immediately 
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imply that the lower bound of the integral is no longer e that we had in the previous 
subsection, but it is 



eur 



e 



(3.41) 



where Umax is the lowest value from the inequality ( p.24| ) that we chose in the previous 
subsection (to avoid clutter we use the same notation). Note that e'is independent 
of ^max- Using this we can now write the distance d between the quark and the 
antiquark as: 



d = 2Voz\ 



5 

max 



dv V 



V^m^ ^ma 
(•^"'^"'^max) 



-1/2 



(3.42) 



where we have defined Vq = ^n-^max ^^'^ sum over repeated index is implied as before. 

Comparing Vq with ( |3.39| ) we see that Vq can be made small if << 1 (which 
is consistent with ( p.9|) and (|3.34| )) as all An for n > 3 are very small. Additionally, 
from (|3.42|) , the term -Zmax^o^ will imply that it will be sufficient to restrict Vq to the 
following series: 



A, 



+ 



A. 



Aa 



(3.43) 



because A^ onwards are very small to consistently maintain the reality of d in ( p.42| ) 
as well as their QsNf dependences from ( p.9|) . This means d in (|3.42|) can be further 
simplified to: 



d 




dv 



^ ^ ^ ^max ^0 



1 - (A - 1^2) 



^max'^ 



32n 



+ 



^ {g, - 4^2) ^^ax +\iG2- 8A2) V^Z, 



2^10 
max 




Since we have taken both A2 as well as B2 to be very small, and plugging in the 
value of Vo from ( p.43| ) we see that d is dominated mostly by inverse Zmax terms, i.e 



d 



3Zn 



0^2 



-OiAl) + O 



(3.45) 



The renormalised Nambu-Goto action on the other hand takes the following form: 
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(3.46) 
1 \ 
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where its clear that the string action is dominated by the inverse ^^max terms. Now 
substituting (|]4|) and ( ^A^ in (U), we get 



QQ 



Arc 



2A2 - Q2 



^4 



Aa 



+ ... 



d + O 



(3.47) 



which is the required hnear potential between the quark and the antiquark. The 
above potential can also be rewritten as: 



Va 



d + O 



(3.48) 



where aru^^ = 



etc. It will soon 



j_ and Ho = = -^,^2 = -^(1 + Al 

become clearer why we want to express the potential ( p. 481) in this way. 

However there is one nagging issue that might be bothering the reader, namely, 
how do we know that the potential ( |3.48|) or equivalently ( |3.47|) only has a linear 
term? To answer this question convincingly, we go to the next section where we 
provide a generic derivation of the linear term. 



3.1.3 Generic argument for confinement 

In the above subsection we argued for the linear potential taking all An for n > 1 to 
be small. This is consistent with the supergravity limit of our background because 
in this limit we expect (^s — ?• and QsNf — > with g^N 00. For these choices of 
An, ( p. 391) will imply Zmax < 1 because we expect Umax to take the highest value in 



Region 3. Under such a situation, condition like ( |3.4CI|) will be fully under control. 



and an analysis of the Wilson loop above reproduces the required linear potential at 
large d. 

However a little thought will tell us that the above derivation cannot be the 
complete story. What if Wmax, in appropriate units, is of order (1 — e) where e — 0? 
In that case its inverse Zmax is of order 1, so both Umax and Zmax can no longer be 
good expansion parameters. We may also consider simultaneously the case where gs 
is no longer small so that An for n > 1 are not small either. Such choices will take us 
away from the supergravity limit that we have been following. In this limit, we want 
to ask whether we can still show linear confinement of quarks. Or more generically 
we want to study confinement for a choice of Wmax that saturates the inequality ( p.28| ) 
but does not presuppose any limiting behavior of Mmax,-^n or Qn- 

In the following therefore we will analyze the integrals ( p.25| ) and ( p.33|) in the 
limit Mmax is close to it's upper bound set by (|3.28|) (see also [0). In particular if 



Umax is the upperbound of Umax, then it is found by solving 

1 



^(m + l)An+3urx' = 1 (3.49) 
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We observe that both the integrals ( p.25 ) and ( 3.33 ) are dominated hj v ^ 1 be- 
haviour of the integrands. Near v = 1 and Wmax — ^ Umax the distance d between the 
quark and the antiquark can be written as: 



d 



A ii" 



dv 



A it" 



VA(l-t;) + B(l-t;)2 
'logA - log (2 v/B(ATB) + 2B + A 



B 



(3.50) 



where note that we have taken the lower limit to 0. This will not change any of our 
conclusion as we would soon see. On the other hand, the renormalised Nambu-Goto 
action for the string now becomes: 

"1 dv ^ 
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where A and B are defined as: 
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B 
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A 



(3.51) 



(3.52) 
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Observe that in the integral (|3.5CI| ) and ( p. 511) we have to take the limit Wmax — ^ Umax- 
So A, B should be evaluated in the same limit. Interestingly, comparing ( 3.52 ) to 
( |3.49|) we see that 



lim 



A ^ 



(3.53) 



thus vanishes when computed exactly at Umax- The other quantity B remains finite 
at that point and in fact behaves as: 
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(3.54) 



Our above computation would mean that the distance d between the quark and the 
antiquark, and the Nambu-Goto action will have the following dominant behavior: 

2v^^mUS^Umax log € 
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(3.55) 
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which means both of them have identical logarithmic divergences. Thus the finite 
quantity is the ratio between the two terms in (|3.55|) . This gives us: 



Cren 
d 



T A ii" 



T X constant 



(3.56) 



Now using the identity ( |3.4| ) and the above relation ( p.56p we get our final result: 

(3.57) 



QQ 



"max 



7ru„ 



d 



which is the required linear potential between the quark and the antiquark. Observe 
that the above potential has exactly the same form as ( p.48| ), justifying the fact that 
0{d'^) terms are not generated for this case. 

Before we end this section one comment is in order. The result for linear con- 
finement only depends on the existence of u^ax which comes from the constraint 
equation ( p.49|) . We have constructed the background such that Umax lies in Region 
3, although a more generic case is essentially doable albeit technically challenging 
without necessarily revealing new physics. For example when u~^j^ is equal to the 
size of the blown up at the IR will require us to consider a Wilson loop that goes 
all the way to Region 1. The analysis remains somewhat identical to what we did 
before except that in Regions 2 and 1 we have to additionally consider i?Ns fields of 
the form and log u respectively. Of course both the metric and the dilaton will 
also have non-trivial w-dpendences in these regions. One good thing however is that 
the Wilson loop computation have no UV or IR divergences whatsoever despite the 
fact that now the analysis is technically more challenging. Our expectation would 
be to get similar linear behavior as ( p.57|) here too. We will however leave a more 
detailed exposition of this for future works. 

3.2 Computing the Nambu-Goto Action: Non-Zero Temperature 

After studying the zero temperture behavior it is now time to discuss the case when 
we switch on a non-zero temperature i.e make g{u) < 1 or equivalently the inverse 
horizon radius, Uh finite in (p. 81), where 



9[U) 



1-^ 



(3.58) 



Choosing the same quark world line ( |3.7| ) and the string embedding ( |3.10D with 
the same boundary condition (|3.11|) but now in Euclidean space with compact time 
direction, the string action at finite temperature can be written as 
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(3.59) 



-37- 



where QmU^ is defined in (|3.19|) and the correction to QmU^ is suppressed by g1 as 
well as u'^/uj^ because the background dilaton and non-zero temperature induces a 
slightly different world-sheet metric than what one would have naively taken. To 
avoid clutter, one can further redefine these corrections as: 



(3.60) 



Minimizing this action gives the equation of motion for u{x) and using the exact 
same procedure as for zero temperature, the corresponding equation for the distance 
between the quarks can be written as: 
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(3.61) 



Once we have d, the renormalized Nambu-Goto action can also be written following 
similar procedure. The result is 
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which is somewhat similar in form with ( p.33|) , which we reproduce in the limit 



Uh ^ CO. Also as in (|3.33|) , we have defined yV^ 

Now just like the zero temperature case, requiring that d be real, sets an upper 
bound to Mmax, that we denote again by Umax, and is found by solving the following 
equation: 



i(m + l)^.,3U-if + ^ n (a: - ^) (%^)' 

J- fc=o ^ / \ h / 



Ami 



2j-l 



: 1 

(3.63) 



Once we fix Uh and the coefficients of the warp factor An, u^ax will be known. We 
will assume that Umax lies in Region 3. 

Rest of the analysis is almost similar to the zero temperature case, although the 
final conclusions would be quite different. To proceed further, let us define certain 
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new variables in the following way: 



l—m ( l—m \ 1 
h \ 4: )■ 



n 

k=0 



2 



I -m> A 



= I - m < A 

= Ai l-m = (3.64) 

As before, we observe that for Wmax — > Umax, both the integrals (|3.61| ),( ^l62D are 
dominated by the behaviour of the integrand near v ^ 1, where we can write 
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where taking the lower limit of the integral to again do not change any of our 
conclusion. On the other hand, the renormalised Nambu-Goto action for the string 
now becomes: 
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where A and B are defined exactly as in (|3.52|) but with An replaced by An given 
by ( 3.64|) above. It is also clear that: 

(3.67) 



lim A 

U m ax ^ U m ax 



and so both ( p.65| ) as well as ( |3.66|) have identical logarithmic divergences. This 
would imply that the finite quantity is the ratio between ( |3.66| ) and ( p.65| ): 
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Now using the identity ( ^.61 ) and the above relation ( ^.56p we get our final result: 
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3.2.1 Analysis of the melting temperature 



To determine the behavior of the potential Vqq as the temperature is raised or 
decreased (or Uh is decreased or increased respectively) we have to carefully analyse 
the behavior of Umax as a function of Uh- Comparing this with ( |3.63| ) and fixed An 
we observe that 



5Ur, 



5uh 



\m{m + l)AmU^ + i n - I) ^Aiu\Aj + /) (| + 2j - l) 

J ■ 111 



-(3.70) 



where the repeated indices are all summed over and the product runs from A; = to 
k = j — 1. Observe that the numerator of ( p.70| ) is always negative, and for large Uh 
the denominator will be positive (because we are taking all An > 0). This means 



(5ui, 



Suh 



< 



and therefore as Uh is decreased (i.e the temperature is increased), Umax increases 

4 

making the ratio ^^^^^ to increase. This in turn would imply that the slope of the 
potential Vqq decreases. Therefore there would be a temperature where the slope 
would be minimum and the system would show the property of melting. 

To start off let us consider (|3.63|) for a simple case where we keep Wmax only to 
quartic order^^. This means: 
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(3.71) 



This is a quartic equation and one can easily solve it for Wmax- To make the analysis 
a little more simpler, let us also assume A3 = 0. Such a choice will immediately give 
us the following potential between the quark and the antiquark: 
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(3.72) 



which tells us that Uh has to be bounded by 

0.84 



Uh > 



AT 



for ( |3.72|) to make sense, and the slope of the potential would decrease as ap- 
proaches this value. On the other hand Wmax increases as Uh is lowered and for 

1.1067 



A 



1/4 



(3.73) 



^■^This is a subtle issue because we are truncating the series ( 3.63 ) and especially to only 
this order to have an analytic control on our calculations. A more generic analysis can be done 
numerically, which we present in the next section. 
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we expect the potental to have a minimum slope where the onset of melting should 
appear. A temperature greater than this is physically not possible because the string 
would break. Also note if we kept 7^ we would have {Ao/AiY^'^ in ( |3.73| ). 

The above conclusion is certainly interesting, but may be a little naive because 
there is no strong reason to terminate the constraint ( |3.63| ) to A4, because ( p.64| ) 
would tell us that higher powers of Wmax will have all the lower ^^'s as coefficients. 
This would make the subsequent analysis complicated, and we have to resort to 
numerical methods. 

On the other hand the above analysis does shed some light on the situation 

4 

where the ratio ^^^^^ > 1. It is clear from our above calculation what happens when 
Uh becomes too small: since by lowering Uh there is an increase in Umax, we are always 
bounded by the constraint: 

Umax < Uh (3.74) 



where the equality would lead to (|3.73|) . Therefore for ( p.72|) and ( p.74| ) to make 



sense, the string connecting quark and the antiquark should break when Umax starts 
exceeding Uh as mentioned above. This is the point where total melting happens, 
and the linear potential goes from the minimum slope a, where 

a = = 1 + 0.32^2 (3.75) 

a 

to zero as soon as the temperature is increased, or alternatively Uh is decreased, 
beyond ( |3.73| ). In fact beyond certain temperature, constraint equation like ( p.63| ) is 
no longer valid, and we are only left with the Coulombic term that eventually dies 
off at large distances. 

3.3 Numerical analysis 

Most of the calculations that we did in the previous sections have been analytic. 
Under some approximations we could see certain important properties of quark- 
antiquark potentials at zero and non-zero temperatures. However as mentioned in 
the footnote earlier, truncating the series ( |3.63| ) as ( p.71D may not capture the full 



story although the above toy example does give us a way to compute the melting 
temperature where the slope of the potential hits a minima ( p. 751 ). Our main conclu- 



sion of the previous section is that there exists a set of warp factors An for which the 
equality Umax = Uh is valid. That means out of a large sets of possible backgrounds 
(classified by the choices of warp factors satisfying EOMs) this equality would se- 
lect a particular subset of backgrounds that allow deconfinement and quarkonium 
meltings at the temperatures ( 3.73|) for a subset of A4 in this approximation. What 



happens if we choose an arbitrary set of warp factors that do not lie in this subset? 
In this section we will perform a numerical analysis to study the behavior of the 
quark- antiquark potentials for this case at all temperatures. 
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For this particular choice of coefficients of the warp factor we can numeri- 
cally compute the interquark distance d in ( p. 251) , (|3.61|) and the NG action S^q 



( p.29| ),( |3l62D , for various values of Umax and using this, plot Vqq as a function of d. 
The analytic behavior of d and Vqq discussed in the previous sections for Wmax very 
large and small will turn out to be consistent with our numerical analysis although 
the property of melting will not be visible now. For simplicity we will choose the 
following values for the coefficients of the warp factor: 

A = 1, ^0 = 1, A2 = = = Qa = 0.24 (3.76) 

with Qs ~ 0.02 and Nf = 24. This is indeed a reasonable choice, and hopefully 
satisfies EOMs despite being outside the required subset, because as we saw from F- 
theory, all corrections to warp factor due to running of the r field comes as 0{ggNj). 
Figure 7 shows how the inter-quark distance d varies with Umax- For T = 0, from the 
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Figure 7: Inter-quark distance d as a function of Umax evaluated for our choice of warp 
factor given earlier. The red curve is the zero temperature limit. Here T = l/n^ = 
henceforth. 



figure we see that there exists an upper bound Umax near which c? — )■ 00. A similar 
analysis also gives as Umax Umax, V'qq — i- 00. By increasing Wmax near Umax, we 
can get all the values of d and Vqq and using this we can plot Vqq as a function of d 
as shown in Figure 8. Note that for large the potential grows exactly linear with 
distance indicating linear confinement. This is consistent with our earlier analytical 
calculations. Additionally for small distances, the potential behaves like a Coulomb 
potential as also predicted by our analysis. On the other hand for T 7^ 0, again for 
the choice of our warp factors ( |3.76|) , from Figure 7 we observe that for every T 7^ 
curve, there exists a maximum value of ci, say rfmax and therefore for every value of c?, 
there are two distinct values for Umax- Such a behavior has also been observed in [pOf 
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Figure 8: Quarkonium potential as a function of inter quark distance d at various tem- 
peratures. Note the linear and the Coulombic behaviors at large and small distances 
respectively. 



for the AdS case and in ||T9| , for the pure Klebanov-Strassler case. This means for 
a particular choice of d and boundary condition u{±d/2) = 0, there are two fZ-shaped 
strings with two values for Wmax, namely Mmax,i and Mmax,2 with 'Umax,2 > Mmax.i- As 
'Wmax,2 > Winax,i, the ?7-shaped string with Mmax,2 has higher energy than the one with 
Wmax.i- We have denoted the [/-shaped string with Mmax.i by branch I and the one 
with Mmax,2 by branch II in Figure 8. It is clear from the plot that branch I has lower 
energy than branch II. Thus at small d, the potential for branch I behaves as Coulomb 
potential and by comparing this to the zero temperature Coulomb behavior, we see 
that the T 7^ Coulomb potential is suppresed. 

Now note that as we lower the temperatue, the value for cimax increases. There- 
fore in the limit T — > 0, cimax 00, which is perfectly consistent with our zero 
temperature curve in Figure 7. This implies as T — > 0, the curves in Figure 7 con- 
verge with the zero temperature curve, which in turn blows up at u^ax- Such a 
behavior is not inconsistent with the high temperature case because we can view the 
T = curve to go straight up and never come down, resulting in a single solution 
for the [/-shaped string for every d. The high temperature curves go upto some dmax 
and then come down. For large d < (imax, we have linear potential for branch I — 
which is suppressed compared to the zero temperature curve as shown in Figure 8. 
As the temperature is increased more is the suppression and smaller is the value for 

In Figure 8 with our choice of the warp factor coefficients the slopes of linear 
potentials and the resulting suppressions are not very significant. For a better view 
of the suppressions, we present a blown up version of Figure 8 in Figure 9. 
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Figure 9: Suppressions at non zero temperatures of the linear potential magnified by 
choosing a slightly different values of the coefficients of the warp factor given earlier. In 
this figure one can clearly see how with high temperatures the quark-antiquark potentials 
melt. 



For d > (imax, there is no real solution to the differential equation for u{x) with 
boundary condition u{±d/2) = and thus there is no [/-shaped string between the 
quarks. Thus for d > (imax, the string breaks and we have deconfined quarks. Also 
as rfmax decreases with increasing temperature, at high temperatures the quarks get 
screened at shorter distances — which is consistent with heavy quarkonium suppres- 
sion in thermal QCD. 

Our numerical analysis and the plots should be instructive for a generic choice 
of the coefficients of the warp factor where we again may not see the melting tem- 
perature. To study the generic case, first consider the zero temperature limit. Note 
that the existence of real positive u^ax is guaranteed if all An in ( p.64| ) are positive. 
However this may not be true if the original warp factor coefficients An are a finite 
set^^. This is because for a finite set of ^„'s, there will be some AnS which are 
negative and the equation 



m + 1 



Am+3'^ 



rn+3 



(3.77) 



may not have any real positive solutions. This would imply the existence of d^^^ 
and consequently two branches of solution. For large d < rfmax? we will have linear 
potentials with suppressions at higher temperatures with lower values of (imax- 

On the other hand, at zero temperature if all An > 0, there is always a positive 
real Umax and we will have linear potential at large distances. If some AnS are 



24, 



As should be obvious from ([3.64|) , a finite set of An still implies an infinite set of An 
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negative, we could have no real positive solution Umax which means < Wmax < oo 
as there is no black hole horizon. In this case the behavior of d will be dominated 
by (i ~ ttmax and that of Vqq will be dominated by Vqq ~ l/'Umax which means 
Vqq ~ — and we will have the Coulomb potential for all d. 

Our above numerical analysis certainly illustrates the decrease in the slope of 
the linear potential as the temperature is increased but does not show us the melting 
temperature. What would happen if we restrict our warp factor choice to the required 
subset of In this case all the high temperature curves will grow linearly and 
will not come back, and at a certain temperature the slope of the curve will drop to 
zero. This will be the melting temperature. In this paper we will not pursue this 
anymore, and more details will be presented elsewhere. 



4. Conclusions and Discussions 

In this paper we have tried to achieve two goals: First is to find the dual to large 
gauge theory that resembles large N QCD i.e. at far IR the theory confines 
and at far UV the theory shows a conformal behavior. We then extend this to 
high temperatures. Second, is to compute the heavy quarkonium potential in this 
theory both at zero and non-zero temperatures. We have shown that, under some 
rather generic conditions, zero temperature linear confinement for heavy quarkonium 
states can be demonstrated. At high temperature, the expected deconfinement and 
quarkonium melting follow from our analysis. 

There are however still a few loose ends that need to be tightened to complete 
the full story. The first one is the issue of supersymmetry. Although we have shown 
that all the unnecessary tachyons can be removed from our picture, this still doesn't 
imply low energy supersymmetry in our model (at zero temperature). Having no 
supersymmetry should be viewed as desirable because we don't expect low energy 
susy in real world! However susy breaking in our model may trigger corrections in 
the potential that need to be worked out. We expect these corrections to change the 
coefficient of the linear term without generating an 0{d'^) term. These corrections 
should be higher orders in QsNj, so will not change any of our conclusions presented 
here. This is because the linear potential arises from the limit where the Nambu-Goto 
action and the distance d exhibit identical logarithmic divergences. This behavior 
should remain the same whether or not we have low enery susy in our model or not. 
Thus the linear confinement argument is particularly robust for our case. On the 
other hand the Coulombic behavior is model independent, so the coefficient of the 
Coulombic term should remain unchanged whether we take susy or non-susy models 
(see for example the model of ||3^). One may choose other embeddings of seven- 
branes, like p2[ or the model studied in |]18[, to study the quarkonium potentials at 



zero and non-zero temperatures. But such choices of embeddings will not change our 
main conclusions. 
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The second one is the issue of Higgsing that breaks the gauge group from SU{N+ 
M) X SU{N + M) to SU{N + M) x SU{N) in the gauge theory side. As mentioned 
in sec. 2.3, the story in the dual gravity side is somewhat clearer. What one needs is 
to analyse the gauge theory operators carefully that will allow the above mentioned 
Higgsing. We leave this for future work. 

The third one is to find the precise set of warp factors An that satisy EOMs 
and allow us to get the melting temperature for the heavy quarkonium states. In the 
previous section we gave a numerical analysis with an arbitrary truncated set of warp 
factors that shows the decrease in the slope of the linear potential with increasing 
temperatures. Our numerical analysis certainly shows the possibility of melting, but 
doesn't tell us the melting temperature. On the other hand our analytic way of 
getting the melting temperature is not very generic. So it would be interesting to 
find the full set of warp factors to complete the story^^. 

Finally, we haven't actually computed the exact gauge fluxes on the seven and 
five-branes that would cancel the tachyons in this model. Following the works of [^| 



this may not look like a difficult task to do, at least for the fiat background. What 
makes it non-trivial here is that all the branes are embedded in a curved background. 
Quantisation of strings in a curved background is highly non-trivial, so it'll be rather 
challenging to work this out in full details. Nevertheless, if we restrict everything 
to Region 3 and away from the brane-antibrane systems, these subtletes will not 
affect our results in any significant way. Happily this is the regime where most of 
our calculations have been performed in this paper. 

Note: As this draft was being written, we became aware of the work of Gaillard et 



al which has some overlap with this paper. See also the earlier work [R3] 
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^Note that in the series of papers 32 similar screening hke what we have in ( 3.69 ) is also 

observed using completely different techniques than ours. Our prediction then would be that the 
square-root suppression that we see in ( 3.69| ) at high temperatures is universal. Of course it would 
be interesting to figure out the Coulomb screening at high temperature also. 
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A. Complete analysis of a background configuration 



In this appendix we will, for illustration, compute the background that appears from 
the backreactions of the seven branes without any three-form fluxes but with five- 
form fluxes. These five-form fluxes are the remnant of the three-branes. So in the 
gauge-theory side we will have a system of D3 branes and Nj seven branes. The 
excitations of these D3/D7 branes are described by a gauge theory with SU{N) x 
SU{N) color symmetry and SU{Nf) x SU{Nf) flavor symmetry. Holography dictates 
that the near horizon geometry sourced by these D3/D7 branes is dual to the gauge 
theory (decoupled from gravity) which arises from the brane excitations. Of course 
there are two ways to get the gravity dual picture from the brane configuration. We 
could first obtain the geometry sourced by the D3/D7 system and then take the near 
horizon limit of it. Or first compute the near horizon geometry of the D3 branes, then 
place seven-branes in that geometry and finally compute the backreaction. Both this 
approaches are identical in the limit where the stack of seven-branes are separated 
from the stack of D3-branes. Therefore for technical simplicity, we will adopt the 
latter approach. We will also put the configuration in a conifold setting to mimic 
what we did in |]13 . 



The near horizon geometry of the stack of D3 branes placed at the tip of a conifold 
is AdS^ X T^'^. We embed the seven-branes in this background with the world- 
volume filling up four Minkowski directions plus a four-cycle in T^'^, and compute 
the backreaction. The supergravity action describing the geometry in Einstein frame 
is 

^suGRA = ^ / d'^'^VG (^R + - ^\Frf^ + S'^!; + j^C^^ A Ri2) 

(A.l) 



where r is the axio-dilaton field, is the five-form fluxes sourced by the D3 branes, 
C4 is the pull-back of the four-form potential, G = ^/det G^,^ with G^,^ being the 
metric, i?(2) is the curvature two-form, and S^j is the local action for the seven- 
branes. There are additional Chern-Simons couplings of the seven-brane gauge fields 
to the background RR forms, but we will ignore them for the time being. They will 
appear later. 

Varying the action with respect to the various fields give rise to the background 
equations of motion. Using our earlier notations, they are given by (|2.27|) (without 
the {p, q) five-brane terms) and ( p.28| ) (again removing the (p, q) five-brane contri- 
bution) for the metric and the five-form respectively. Our metric ansatze remains 

1 



ds^ = —= 
Vh 



- g{r)dt^ + dx^ + dy"^ + dz^ + Vh g{r) ^g^dr'^ + gmndx'^dx^ (A. 2) 



as before, with g{r) being the Black-Hole factor and h being the warp factor that 
depends on all the internal coordinates (r, ^j, 0j, ip). To zeroth order in ggNj we have 
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our usual relations: 



= ^, g^'^ = 1 - 4> 9f} = 1, gtidxV = dsU (A.3) 



^4 - lyH 



But in higher order in QgNf, both the warp factor and the internal metric get modified 
because of the back-reactions from the seven-branes. We can write this as: 

h = /.M + /.w, g„ = g^^} + gl"}, gmn = g^l + gl^n (A.4) 

where the superscripts denote the order of gsNf. 

Using the full F-theory completion of the background as discussed in section 2.3, 
we know that near any one of the seven-branes, i.e z ^ Zk'. 

T ~ \og{z-Zk) (A.5) 

and therefore the internal metric components will typically behave as: 

V ^[1] [1] ^ V- [1] log^(r) 



with a[)^„ • • independent of r but depend on the internal coordinates (■?/', 6i 



Now away from z ~ z^, we can use similar discussion as in sec. 2.3 with axio- 
dilaton r behaving exactly as ( p.33| ). In that case the internal components of the 
metric become: 



oo [1] oo [1] 

V \ ^rr^V^ ^ ril V 

i=0 j=0 



where again a|^'„ ^ are independent of r but depend on the internal coordinates. 

Now to find the precise expression for /i, we use the five-form equation ( |2.28| ). 
Since both the three-forms vanish for our case, and additionally if we embed part of 
the spin connection in the seven-brane gauge connection with appropriate number 
of antiseven-branes, we can easily derive the following warp factor equation: 

E U-g''^^^.^^ = sources (A.8) 

as before; where (m, n) now run over all the internal six coordinates. The generic 
solutions of all the above equations would be 



^4 



J.4 



i=l 



for large r 

for small r (A. 9) 



- 48 - 



which is of course the solutions (|2.7|) discussed before, and away from the interpo- 
lating region, i.e Region 2. In the next appendix we will discuss a way to determine 
the coefficients and bij. Note that the above form of the large r warp factor also 
implies that the effective number of colors is given by: 



Thus keeps growing with radial coordinate r for all < and A^^efr = iV at the 
bounday r = oo. Alternatively, with a change of coordinate u = 1/r, the metric will 
take the following form: 



which is the metric ( p.8D and dsj^^ is the metric of the deformed T^'^. We note as 
before that the coefficients C„ depend on the coordinates {tp,6i,(f)i), the locations 
of the seven-branes, and the number of colors and flavors, and Nf respectively. 
For all Cn > 0, Ncs grows with decreasing u with maximum value at the boundary 
u = 0. 



B. Solution to Einstein equations 



One last thing that needs to be studied is the deviation of the internal metric from 
the usual Ricci-flat metric. We already hinted this in sec. 2.3, but now we will keep 
all the angular coordinates and see whether we can flnd a possible solution. Our aim 
therefore is to solve (|2.27|) , keeping only the axio-dilaton sources ( p. 331) and not the 
seven-branes and antiseven-branes local energy- momentum tensors. A justiflcation 
for this can be easily provided: we are using the efective embeddings (|2.8|) and ( p.37| ) 
for branes and antibranes respectively. The axio-dilaton fleld ( p.33|) can alternatively 
be written as: 



r 



cosec 



cosec 



(B.l) 



The above behavior, as discussed earlier, is valid for all r (in Region 3). Therefore r 
has a nice Taylor series expansion and the right hand side of ( |2.27|) will also have a 
smooth Taylor series. As an illustration, consider 



dr dr 



dr dr ^ 



k,l 



r 



*i+02-^)(fc-Ocosec^+' 



— I cosec 



2 



,^^2)(B.2) 
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mny^i, 62) depends on the partial derivative \drT\ ■ As the source is a Taylor 
^ we must have Ricci-tensor Rrr to be a Taylor series. This in turn would 



where F, 
series in r 

imply that the unwarped metric components should go like: 



9r^ 



E 



(B.3) 



where j = i — 2 or j = i depending on whether we are choosing the radial direction 
or not (see sec. 2.3 for details). The coefficients amni = o-mnii^, 4>^4') i-6 functions of 
all the internal coordinates. 

Now to explicitely solve for the partial differential equations involving amni we ob- 
serve that the right hand side of ( p. 27] ) as a function of internal coordinates {ip^ipi^Oi) 
is of the form e^i<l>i+<l>'i-'^) and cosec (^). So the most general form for amni should 
also involve these variables^^. Thus we introduce variables: 



Z2 



COS 



COS 



— + 2sm — 
2 / V 2 



+ zsin 



(B.4) 



and using which we make the following ansatz for 



+ 



k,l,k+l=a 



~i]op 



+ 



+ 



ZiZ j ZqZ^ 



Zii^i^. 



■Zi2a~-lZi2a 



(B.5) 



where summation over repeated index is assumed. Note that with this ansatz written 
in terms of coordinate Zj, the differential equations in (|2.27|) become a set of algebraic 
equations because sec'^(6'j), cosec"(6'j) etc. can be written as a linear combinations 
of ij. Then it is straight forward to equate coefficients of various powers of Zi and 
obtain the constant coefficeints aiij2..i2n- This gives us the precise procedure to obtain 
the exact solutions for all amni- In Region 2 we have to be more careful because there 
are, in addition to the axio-dilaton, other sources. But since all our calculations are 
restricted to Region 3, we need not worry too much about the non-Ricci flatness of 
the internal manifold. 
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